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f H {x ; d) = 2 hmsup -^ . 

tiO.d'-xi I 

wj Riemann's second upper derivative of f at x a in the directions d,w G R n 
is defined as: 

-a f(x + td + tw) - f(x + td) - f(x + tw) + f(x ) 
f R (x ;d,w) = hmsup -= . 

no t 

iv) Dini's second upper derivative of f at xq in the directions d, w G R™ is 
defined as: 

-„ Vf(x + td)w-Vf(x)w 

f d (xq; a, w) = limsup . 

tj.o t 

v) Yang-Jeyakumar's second upper derivative of f at x a in the directions 
d,w G R" is defined as: 

-fii , ,. V/(#o + td+ tz)w — V/(#o +tz)w 
fy\ x o\ u , v ) — SU P hmsup . 

2£R" tj.0 t 

vi) Clarke 's second upper derivative of f at xq in the directions d,w G R" is 
defined as: 

-j, X7f(x' + td)w-X7f(x')w 
} c (x ;d,w)= hmsup . 

x'->x, ij.0 t 
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vii) BenTal-Zowe's second upper derivative of f at xq in the directions d,ic£ 
K™ is defined as: 

-,, f(x + td+2- 1 t 2 w)-f(x )-t\7f(x )d 
f B (x ;d,w) = 21imsup - . 

Remark 2.1. In an analogous way one can define lower derivatives and we 
denote them by f . 

In the following results we will give some relations among the derivatives of 
definition 2.2. 

Remark 2.2. In [4-9] is given the following chain of inequalities: 

J D (x ; d, w) < J Y (xo; d, w) < J c (x ; d, w). 

Furthermore f Y {xo;d,w) = f c (xo;d,w) if and only if the map f Y {-',d, w) is 
upper semicontinuous . 

Remark 2.3. In [9] is given the following characterization of Clarke's gener- 
alized derivative: 

—r-d.W . . 

f c (x ;d,w)= hmsup 

where: 

Af w f(x; s,t) = f(x + sd + tw) - f(x + sd) - f(x + tw) + f(x). 

From this characterization one can trivially deduce that f b (xq', d, w) < fci x o', d, w). 
Proposition 2.1. Let f be a function of class C > at xq. Then: 

J B (x ;d,w) = Wf(x )w + J P (x ;d). 
Proof. Let l B e d B f(x ;d, w) and l P e dpf(x ; d) where: 
,2^_.^^_/,_o lim /(go + *g* + 2 ~ H l w ) - f( x o) - tkVf{x )d 

k — *+oo tu 



d B f(x a ;d,w)=^l = 2 : lim ^-^ ^ * , 2 , *fc I 

t k 

& f( « J, o v f{xo + t k d) - f(x ) - t k Vf{x )d 
d P f{x ; d) = <l = 2 lim -, ,t k J. 

[_ fe— >+oo t k 

Eventually by extracting subsequences, we have: 

, _ .. f(x + t k d + 2- 1 tlw)-f(x )-t k Wf(x )d 
l B = 2 lim -2 = 

k — )-+oo tu 

/(.to + tkd + 2- 1 t 2 k w) - /(xq + t k d) t f(x + tkd) - f(x ) - t k Vf(x )d 
Z lim 9 -\- 9 

fe^+oo t% t% 



f(xp + t k d) - f(x ) - t k Vf(x )d 

fe^+oo ' t. 



2 ? hm V/((t)wH -J = V/(a;o)w + /p 



D 
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Lemma 2.1. Let f be a function of class C > at Xq. Then f P (xo',d) — 
J H (x ;d) and f' p {x \d) = f' H {xa\d). 

Proof. It is trivial that f P (xo;d) < f H (xo;d). Vice versa, let I e djjf(xo;d) 
where: 



J(x + tkdk) - f(xo) - t k Vf(x )d k 

t 2 

'■k 

Then there exist t k J. 0, d k — > d s.t.: 



d 2 H f(x ; d) := {/ = lim 2^^ R/ J y J y v ' K - ,t k | 0,d k -> d}. 



' r f(x + t k d k ) - f(x ) -t k \7f(x )d k 
- = urn o = 

2 fc — >+oo ^i 

,. /(zo + ifcdfc) - /(^o + t k d) + f(x + t k d) - f(x ) - t k Vf(x )d - t k Vf(x )(d k - d) 

lim „- . 

k — >+oo t. 

Taking eventually a subsequence: 



f(x + t k d) - f(x ) - t k s/f{x )d _V , -ft . 
72 — y ' ^ Jp\ 2 

l k z 

f(x + t fc rffc) - /(.tq + tfecQ - t k Vf(x )(d k - d) 
t 2 

'k 

tfcV/(&)(d fc -d)- t k Wf(x )(d k - d) 



lim 

k — >-+oo 



lim 

k — *+oo 



< 



liu) K\\d k - d\\\\Z k - x \\ =Q 



k — >+oo t k 

where £ k e [a?o + tkdk,xo + t k d] and then t k ~ x ° — » d. Therefore: 

,. f(x + t k d k ) - f(x + t k d) -tkVf(x )(dk-d) 
lim s = 

fc^ + OO tf, 

and then V = I. Then f h(xq', d) < f p (xq', d). The proof of the second equality 
is analogous. □ 

Theorem 2.1. Lei / 6e a function of class C 1 ' 1 ai xo- Then: 

i) 7ffOo;d) =Jp{xa;d) < J D (x ;d,d) < J Y (xo;d,d) <J c (x ;d,d). 
ti) Jh( x o; d) = J P (x a ;d) < J R (x ;d, d) < J Y ( x o; d, d) < J c ( x o; d, d). 

Proof. i) From the previous remarks, it is only necessary to prove the in- 
equality / p(xq', d) < f u(xo', d, d). If we take the functions 4>i(t) = f(xg + 
td) — iV f(xo)d and <p2(t) = t 2 , applying Cauchy's theorem, we obtain: 

f(x + td)-f(x )-tvf(x Q )d = Mt)-MQ) = 
' ~ ~ t 2 "~ ' Mt)-M0) 

MZ) ^ Vf(x + Sd)d-Vf(x )d 
where £ = £(£) e (0, t), and then f P (xo\ d) < frj( x o] d, d). 
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ii) From the previous remarks, it is only necessary to prove the inequalities 
f P (xo;d) < f R (xo;d,d) < f Y (xo;d,d). Concerning the first inequality 
from the definition of / P (xo; d) we have: 

t 2 — 
f(x + td) = f(x ) + tVf(x a )d + -f P (x ;d) + g(t) 

where: 

lim sup — =- = 

no * 

and: 

f(x + 2td) = f(x ) + 2tVf{x )d + 2t 2 J P (x ;d) + g(2t) 



Then: 



where: 

g(2t) g(2t) 
limsup — 7T— = 4hmsup =- = 0. 

no t2 no 4t 2 

f(x + 2td)-2f(x + td) + f(x ) 
lim sup ^ = 

no t z 

r t 2 f P {x ;d) + g{2t)-g{t) 

lim sup = > 

tio * 

-i" , a ,. 5(2t) ,. 5W 

f p(xo; d) + limsup '■ — limsup — — . 

no * no * 

Then f r (xq; d, d) > f p (xq; d). For the second inequality we define 4>\{t) = 
f(xo + 2td) — 2f(xo + td) and <f>%{t) = t 2 . Then, by Cauchy's theorem, we 
obtain: 

/(.to + 2td) - 2/(2:0 + td) + /(gg) = 0i(t) - 0i (0) ^ 

t 2 M*)-Mo) 

0j(Q _ V/(x + 2grf)rf - V/(a;o + tfrflrf 

where £ = £(£) € (0, t), and then fp(xo; d, d) < f Y (xo; d, d). 

□ 



3 Optimality conditions for constrained optimiza- 
tion problems 

The aim of this section is to give necessary and sufficient conditions for the 
following optimization problem: 

OP) min/(a;) 

where S is the feasible region defined as: 

^{leM™: hi{x) =0,i = l...m,gj{x) < Q,j = l...r,x€ X}. 
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We suppose that the functions / : E™ -> E, h t : E™ -► E and 5j : E" -> E, 
are of class C 1,1 , i = 1 . . . m, j = 1 . . . r and X is a closed convex subset of E™. 
From these conditions one can easily deduce the results in [19, 28, 49]. 

The following definitions will be useful in order to obtain optimality condi- 
tions for OP). 

Definition 3.1. Let A be a subset o/E™ and xo £ clA, where clA denotes the 
closure of A. The sets: 

• i) WF(A, x ) = {d : 3t k J. 0, x + t k d £ A, Vk} 

• %%) F(A, x ) = {d : Vi fe J. 0, x a + t k d £ A, Vfc} 

• Hi) T(A, x ) = {d : 3d k — ► d, 3t k J. 0, x + t k d k £ A} 

• iv) T 2 (A, x , d) = {w £ E n : 3w k -► w, 3i fe J. 0, x + i fe rf + 2- Y t\w k £ A} 

• v) T 2 (A,x a ,d) = {w £ E" : 3 Wfc -» w,3t fc I 0,3 7fc I 0,7^ _► 0,a; + 

are called, respectively, the cone of weakly feasible directions of A at Xo, the 
cone of feasible directions of A at x , the contingent cone of A at Xq, the second 
order contingent set and the asymptotic second order contingent set of A at xo 
in the direction d. 

Theorem 3.1 (Proposition 3.3.11, [3]). Let xq be a local minimum point of 
the problem OP) where f, hi, and gj are of class C 1 from E" to E and X is 
a closed convex set. Then there exists a scalar /j,q and Lagrange multipliers \ 
and /ij satisfying the following conditions: 

• i) for all x £ X , we have: 

* m r \ 

MoV/(.t ) + ^2 K^hi(xo) + ^2 Hj^9j( x o) ix - xo) > 
»=i j'=i / 

• ii) A*j > , j = . . . r 

• Hi) Ho, A, and jjlj are not all equal to 0. 

• iv) in every neighborhood N of xo there is an x £ X ("1 N such that 
\ihi(x) > for all i with Xi ^ and fijgj(x) > for all j with fij =/= 
(this condition implies complementary slackness (ijgj(xo) = 0-) 

Theorem 3.2 (Proposition 3.3.12, [3]). Let xo be a local minimum of the 
problem OP) and assume that the following two conditions hold: 

• a) there does not exist a nonzero vector A = (Ai, . . . , A m ) such that: 

( Yl A ' Vfe i( s o) ) ( x - x o) >0,Vi€l 
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• b) there exists a feasible direction d of X at xq such that: 

Vhi(xo)d = 0, % = 1 . . . m, Wgj(xo)d < 0, Vj G ^4(a;o) 

where A(xo) = {j ■ Sj( x a) = 0}- 

TTien i/ie previous condition holds with fiQ = 1. 

Given A € M m , /io G IR and [i G R r , let L(x,X, /i) the usual Lagrangian 
function, that is: 

m r 

»=i j=i 

and G(fi) — {x G X : fj,g(x) = 0, h(x) = 0}. We will write L(x) instead of 
L(x,X,fi) when this doesn't create confusion. The function L(-,A,/i) is of class 
C 1 - 1 . 

Theorem 3.3 (Necessary optimality condition). Lei xo &e a local mini- 
mum point and assume that the previous two conditions a) and b) are satisfied. 
Then: 

• there exists A G W 71 , /i G R r , fij > 0, j = 1 . . .r, such that /j,jgj(xo) = 
and V x L(x 0l A, [i)d > 0, Vd € T(G{h),Xq); 

• Vd G T(G(fi), x ) such that V x L(xo, A, ^t)(i = 0, we have: 

V x L(x , A, ^)w; + L xH (x , A, /u; d) > 0, 
VweT 2 (G(^),a;o,d). 
Proof. • i) See theorem 3.1. 

• ii) Suppose, ab absurdo, 3d G T(G{(j),Xo) with V x L(a:o, X,n)d = 0, and 

w; G T 2 (G(\),x a ,d) such that: 

V x £(£n, A, (j,)w + L xH (x , A, /u; d) < 0. 

Since ty G T 2 (G(n),xo,d) then there exist w k —* W, ife j such that 
xq + ifcfi + 2~ 1 t k w k G G(fi). So, by extracting eventually subsequences, 
we have: 

, r £Oo + ifc(d+2~ : 4 fe Wfe)) - L(x ) -t k V x L(x )(d+2- 1 t k w k ) 
I = lim 2 2 

fc^+oo i, 

where I G djjL(xo;d). We observe that: 

L(a; + t fe (d + 2 _1 t fe w fc )) - L(xo) = 

f(x Q +t k (d+2~ 1 t k w k ))+iig(xo+t k (d+2~ 1 t k w k ))+\h(xo+t k (d+2~ 1 t k Wk))- 

f(x ) - ng{xa) - Xh(x ) > 0. 

Then: 

< L(x + tk(d+2- 1 t k Wk)) - L(x ) = 
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t k W x L(x )(d + 2 1 t k w k )+t 2 k -+o{tl) = 

2~ 1 4(V x L(x )w k + I) + o(t 2 k ). 
Taking the limit when n — > +00, we have: 

< \7 x L(x )w + I < V x L(x Q )w + L xH (x ; d). 



□ 



Corollary 3.1. Let x a be a local minimum point for the problem: 

min f(x). 

xGX 

Then: 

i) V/(x )d > 0, Vd e T(X,x ). 

n) Jh(x ; d) > 0, Vd G WF(X; x ), Vf(x )d = 0. 

Proof, i) Trivial, ii) If d G WF(X, x ), Vf(x )d = 0, then G T 2 (X, x , d) and 
soJ H (x ;d) > 0. D 

Theorem 3.4 (Sufficient optimality condition). Let xq E S, fi E R r , 
Mi > 0; IJ>i9i{%o) — 0, A G M m , suc/i i/iai V x L(a;o, \,n)d > 0, Vd G T(S,xq). 
Suppose that Vd G T'(5', xo) smc/i i/ioi V x L(xo, A, /i)d = 0, we have: 

V x L(x )w + L x " h (xq, A, /u; d) > 

Vu> G T (S, xo,d), wd = and: 

V x L(x , A,/i)w > 

\/w G Tq(S, Xo,d), to/ 0, wd = 0. TTien xq is a strict local minimum point. 

Proof. Ab absurdo, let x k — ► Xrjj x^ G 5, be a sequence such that /(xfc) < /(xo), 
x k 7^ xo- If ife = ||xfc — xrj|| and d k = £jT (x& — xrj) then, eventually by extracting 
subsequences, t k I 0, d k -^ d and d G T(S,Xo)- So: 

> f(x k ) + ng{xk) + Xh(xk) - f(xo) - Mfl'(^o) - Xh(x ) = L(x k ) - L(x ) 

and then V x L(xo)d < 0. So V x L(xo)d = 0. Now let s k = \\d k — d\\ and 
w k = sj~ (dfe — d). If dfe 7^ d, Vn G N (eventually by extracting a subsequence) 
then x k = xq + ifed + t k s k w k , s k 7^ and tot ->id / 0. Since d kl d E S 1 , then: 

l|dfe|| 2 = ||d|| 2 + ||s fe w fe || 2 + 2s k w k d 

and then w^d = — 2~ 1 Sfc||u>fc|| 2 . Taking the limit when fc — > +00, we have 
wd = 0. So, by C 1,1 regularity, we have: 

L(x + tfedfe) - L(xo) - t k V x L(x )d k 
I = 2 lim s 
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and then: 



and so: 



> L(x k ) - L(x ) = t k V x L(x )(d + s k w k ) + t 2 k - + o(t\), 



> 2s k t k 1 V x L(x )w k + I + e fe . 
Now if t k s k — > r > then, taking the limit when n — > +oo, we obtain: 

> 2rVi(a;o)u' + Lj' H (x \ d) 

which contradicts the hypothesis, since 2rw G T 2 (S,Xo,d). If tl s k — ► +oo 
then w G Tq(5, xo,cf) and: 

> Vi(x )wA; + tks^ 1 - + ifeSfe 1 o(tfe) 

and so ^ x L(xq)w < 0, which contradicts the hypothesis. If dfc = rf, Vn G N, 
then d £ F(S, Xq) and then € T 2 (S, Xo,d). Following the same calculations of 
the previous case, we obtain a contradiction. □ 

4 Final remarks on optimality conditions for C 1 ' 1 
vector optimization problems by scalarization 

In this section we briefly discuss how to obtain optimality conditions for C ' 
vector optimization problems by scalarization techniques. Many different ways 
to reduce a multi-objective problem to a single-objective ones have been pro- 
posed; we focus the attention on linear scalarization methods which preserve 
the C 1 ' 1 regularity of the involved data. A vector function / : W 1 — * W" 1 is said 
to be of class C 1 ' 1 if and only if each component of / is of class C 11 . Given 
/ : R n -*■ R s , g : R n -» R r and h : R n -» M m , and a given subset I of R", 
consider the following optimization problems: 

1) min/(a;) 

2) min 8f(x) 

where S, the feasible region, is defined as: 

^{leM™: h t (x) =0,i = l...m,gj(x) < 0, j = l...r,xe X}. 

We remember the following definition of vector minimum point and weak vector 
minimum point of the problem 1). 

Definition 4.1. A point x G S is said to be a vector minimum point of f if 
there is no x G S , x ^ x , such that f{x) — f(x) G K™. 

Definition 4.2. A point x G S is said to be a weak vector minimum point of f 
if there is no x G S , such that f{x) — f{x) G intM 111 . 

The following theorem gives a necessary and sufficient condition for vector 
minimum points. 
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Theorem 4.1. [33] Let 9 € intIR™. Then, x is a vector minimum point of 1) 
if an only if it is a minimum point of the single objective problem: 

3) min 6f(x) 

where S' = {x G S : f(x) - f(x) G R™}. 

A sufficient condition for vector weak minimum point is given in the following 
result. 

Definition 4.3. [33] The function f : R n — ► W m is said to be convexlike on 
S C R n if and only if for any X\,Xi G S and any t G [0, 1] there exists X3 G S 
such that: 

tf( Xl ) + (1 - t)f(x 2 ) - f(x 3 ) G R™ 

Theorem 4.2. [33] 

• Let 9 G K™ and x be a minimum point of 2). Lf 9 ^ 0, then x is a vector 
weak minimum point of 1). 

• Suppose that f is convexlike on S. Lf x G S is a vector weak minimum 
point of 1), then there exists 9 G M™\{0} such that x is a minimum point 
of 2). 

Nonlinear scalarizations of vector problems don't preserve the C 1 ' 1 regular- 
ity, that is starting from a C 1 ' 1 vector function the scalarized vector function is 
not a C ' function. Indeed, by linear scalarization, one can use the results of 
the previous sections for the problem 2) in order to obtain optimality conditions 
for the problem 1). 
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Abstract 

In this paper, we investigate the Hyers-Ulam stability problem for the difference equation 
f(x + p,y + q) + ip(x, y)f{x, y) + ip(x, y) = and its applications. 
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1 Introduction 

In 1940, S. M. Ulam [16] raised a question concerning the stability of group homomorphisms: 

Let G\ be a group and let G2 be a metric group with the metric d(-,-). Given 
e > 0, does there exist a 5 > such that if a function h : G\ — > G<i satisfies 
the inequality d{h{xy),h{x)h{y)) < 5 for all x, y € G±, then there exists a 
homomorphism H : G\ — > G2 with d{h{x) , H \x)) < e for all x G G\l 

In other words, we are looking for situations when the homomorphisms are stable, i.e., if 
a mapping is almost a homomorphism, then there exists a true homomorphism near it. If the 
answer is affirmative, we would call the equation of homomorphism H{xy) = H(x)H(y) stable. 
During the last decades, the stability problems of several functional equations have been extensively 
investigated by a number of authors [2-12, 14, 15]. 

Mathematical computations frequently are based on equations, called difference equations or 
recurrence equations that allows us to compute the value of a function recursively from a given set 
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of values. These equations occur in numerous settings and forms, both in mathematics itself and in 
its applications to statistics, computing, electrical circuit analysis, dynamical systems, economics, 
biology, and other fields. The study of linear difference equations is important for a number of 
reactions. Many types of problems are naturally formulated as linear equations [1, 13]. 
In the next section we consider the following difference equation of the form 

f(x + p,y + q) + <p(x, y)f(x, y) + ip(x, y) = 0, (1.1) 

f(x + p,y + q) + ip(x,y)f(x,y) = (1.2) 

with restricted conditions of <p(x,y) and ifj(x,y). 

It is important to provide methods and suitable criterion that describe the nature and behavior 
of solutions of difference systems, without actually constructing or approximating them. In contrast 
with differential equations, since the existence and uniqueness of solutions of discrete initial value 
problems is already guaranteed, one of the problems is the study of asymptotic behavior of solutions 
of the difference system [1]. 

Apart from the stability of solutions of the difference equation, in this paper we examine the 
situations that the difference equation (1.1) is stable in the sense of Hyers and Ulam, i.e., for a given 
(5-approximate function we construct a true solution of the difference equation near it. Throughout 
this paper, let 5 > and p, q G N be fixed, and N denote the set of all positive integers and for 
some nonnegative integer k, N& := {k, k + 1, k + 2, ■ ■ ■ }. 

2 Main Results 

Before taking up the main subject we point out the following situation which is similar to that of 
elementary homogeneous linear differential equation. That is, if a particular solution f p of (1.1) is 
given, then the general solution / of (1.1) has the form f = fh + fp, where fh is a solution of (1.2). 
In the next theorem, let two functions ip : N& x N& — ► (0, oo), ip : N& x Nk — ► K satisfy 

oo j -. 

e{x,y) := VTT , —. — tt < oo, (2.1) 

^ f = o f(x + W,y + iq) 

£ {x,V) ■= 7^77 i ; : tt < °° ( 2 - 2 ) 

j=o \Yi=oVi.x + ip,y + iq) 

for all x,y £ N^. We now investigate the Hyers-Ulam stability problem for the equation (1.1). 
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That is, the difference equation (1.1) is stable in the sense of Hyers and Ulam under the conditions 
subject to (2.1), (2.2). 

Theorem 2.1 Suppose that functions f, ip : J% xNj-»I and (p satisfy the inequality 

\f(x + p,y + q) + ip(x,y)f(x,y) +ip(x,y)\ < 5 (2.3) 

for all x,y £ Nk- Then there exist unique functions T,Th,T p : f% x Nk — ► K swc/i £/m£ T, T p satisfy 
the equation (1.1), T^ satisfies the equation (1.2) and the relations 

\f(x,y)-T(x,y)\<Se(x,y), (2.4) 

\f(x,y) -T h (x,y)\ < 8e{x,y) + \e'(x,y)\, 

\T p (x,y)\ < \e'(x,y)\, 

T(x, y) = T h (x, y) + T p (x, y) 

hold for all x,y £ N&. 7/liminf{/(x,y)|x,2/ 6 N) > 0, then the range of T^ is in (0,oo). 

Proof. We prove the theorem by dividing into various steps. 
Step 1. We show that there exists a limit function 

f(x + 2mp,y + 2mq) ^ • ip(x + jp,y + jq) 
T 2m {x,y) := 2m-1 : — - > (- 1 ) ^n ; : TT 

lli=o v>(& + »p, y + *?) j=0 lIi=o v( x + *P> y + *«) 

for any x, j/ £ J%. 

Replacing x, y by x + p, y + q, respectively, in (2.3), we have 

\f(x + 2p,y + 2q) + ip(x + p,y + q)f(x + p,y + q) + ip(x + p,y + q)\ < 5 (2.5) 

for all x, y £ N&. Combining the last inequality with (2.3), we get the relation 

\f(x + 2p,y + 2q) - ip(x + p,y + q)<p(x,y)f(x,y) + ip(x + p,y + q) (2.6) 

- <fi(x + p,y + q)ip(x,y)\ < 5 + 5cp(x+p,y + q) 
for all x, y £ N&. Now we use the induction on m £ N to prove 

2m- 1 

f(x + 2mp,y + 2mq) - ]~[ cp(x + ip,y + iq)f(x,y) (2.7) 

i=0 
2m— 1 2m— 1 



^ (-i) J V(^ + j>,y + j'9) {J ¥>0* ; + *P)y + 



2m- 1 2m- 1 



< (5 Y^ II p( x + i P>y + i( i)i 

3=0 i=j+l 
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and 

2m 

f{x + (2m + l)p, y + (2m + l)g) + [J <^(x + ip, y + iq)f(x, y) (2.8) 

j=0 

2m 2m 

+ ^2(-l) J ip(x + jp,y + jq) Yl ¥>(x + ip,y + iq) 

3=0 i=j+l 

2m 2m 

<<^xi n ^ + *p> y + *«) 

for all x,y S N/%, where IX?(-) = ^ conveniently if z > j. The validity of (2.7) and (2.8) is easily 
proved by joining it to the other together with (2.3). Thus we obtain the inequalities which play 
an important role in the sequel, 



2m- 1 



f(x + 2mp,y + 2mq) ^ ■ i/j(x + jp,y + jq) 

f{x,y) - 2^ (-1) 



Ui=o v(x + ip,y + iq) 

2m- 1 



j=0 



nLo <p( x + ip,y + k) 



< 



*E 



1 



o Ili=o ^( x + *P> 2/ + *9) ' 



and 



2m 



/(x + (2m + l)p,y + (2m + l)g) y^ ■ 4>(x + jp,y + jq) 

X\f=Qtp{x + ip,y + iq) ' =0 Y[ 3 i=Q v(x + ip,y + ic 



2r« 



7= o Ili=o ^( x + *p. y + *«) 



We claim that the even sequence 



/(x + 2mp, y + 2mg) 



2m- 1 



T2m{x,y) - TT 2,,,-;i . / „ ■ --,-,., , . . . 

lli=o <p{x + ip,y + iq) jr^ Ui=o'P(x + ip,y + tq) 






i>{x + jp, y + jq) 



is a Cauchy sequence. Indeed, for m > n we get from (2.9) 

i 

|r 2m (x,y) -T 2n (x,y)\ 



f(x + 2mp, y + 2mq) 



Ui=o p(x + ip,y + iq) Ui=2n <-p(x + ip,y + iq) 
-/(:,: + 2„p, V + 2n«) - eVi)' -^ V"'^ ~ ; 

j=2n Ui=2n fi X + W,V + W 

c- 2m-l 



< 



E 



IL=o f( x + ip,y + iq) j=2n l\i=2n <P(x + ip,V + iq) 
as n — ► oo. 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



Therefore, we can now define a function T : N& x 1% — ► M. by 



T(x, y) = lim T 2m (x, y) 



(2.13) 
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for any x, y G f%. 

Step 2. We claim that T satisfies the equation (1.1). 
We denote the odd sequence in (2.10) by 
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S 2 m+i(x,y) 



f(x + (2m + l)p, j/ + (2m + l)q) 



2m 



+ »-*) 



^(x + jp, y + jg) 



(2.14) 



Il£i<p(x + ip,V + iq) ^ ' U.i =0 ip(x + ip,y + iq) 

Replacing x, y by x+2mp, y + 2mq, respectively, in (2.3) and dividing by Yli=o V 9 ( x + Wi V + *9)> 
we obtain 

<5 



|<SWi(x,y) + T 2m (x,?/)| < — 2^- 

ll i=0 93(x + ?p,y + zg) 



which implies that a function 5 given by 



S(x,y)= lim S2m+i(x,2/) 



(2.15) 



(2.16) 



is defined since the right hand side in (2.15) tends to as m — > oo, and so S(x, y) + T(x, y) = for 
any x,y G N fc by (2.15). 

Since T 2m (x+p,y + q) = (p(x,y)S 2m+ i(x,y) -ip(x,y), we have 



T(x +p,y + q) = <p(x, y)S{x, y) - i/)(x, y) = -<p(x, y)T(x, y) - ^(x, y) 

for any x,y G Nfc. That is, T is a solution of the equation (1.1). 
Step 3. We conclude that there exists a limit 

f(x + 2mp,y + 2mq) j 



(2.17) 



T h (x,y) := lim 



iii=o 1 <p( x +w,y+ n) 



which satisfies the equation (1.2) for any x, y G f%. 
From (2.9) we obtain that for any x, y G N& 
/(x + 2mp, y + 2mq) 



Ui=o x <p( x + w,y + k) 



f(x,y) 



(2.18) 



2m- 1 



s*E 



1 



j=o n<=o vte + *p> y + *«) 



+ 



2m- 1 

E(-d 

i=o 



ip(x + jp, J/ + jq) 



nLo ^ + ^> y + *«) 



Also from (2.10) we have that for any x,y G N& 

/(x + (2m + l)p, y + (2m + l)q) 



2m 1 

<<5V— - + 

j= o Ui=o <f( x + *P,V + %q) 



+ f(x,y) 



(2.19) 



2m 



E 



(-ipV(^ + j>,y + j'9) 



o nLo vte + iPiy + n) 



(2.20) 
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Using the similar argument to that of (2.12), we obtain that the sequence 

{f(x + 2mp, y + 2mq) 1 
Ili^cr 1 <P(.z + ip,y + iq) J 

is a Cauchy sequence due to (2.18), and the function 7), given by 

/(a; + 2mp, y + 2mg) 

is defined for any x, y € N&. 

As in the case of (2.15), we have from (2.3) 

f(x + (2m + l)p, y + (2m + l)q) /(a + 2mp, y + 2mg) 22 

Wl=Q^{x + ip,y + iq) \^=q X ^p{x + ip,y + iq) 



T h (x, y) = Jim < ^_T r:^ 7\ !> ( 2 - 21 ) 



■0(x + 2mp, y + 2mq) 



< i 

2m 



IL^o <p(x + w,y + n) Wi=o <p(x + ip,y + iq) 

which yields that a function Sh given by 

m _oo y X\l=^{x + ip,y + iq) J 

is defined and Sh(x, y) + T/^x, y) = by (2.22), and hence 

7fc(x + p, y + <?) = v?(x, y)S h (x, y) = -^(x, j/)T h (x, y) 

for any x, y £ N&. That is, T^ is a solution of the equation (1.2). 

Step 4. We obtain that T p (x,y) := T(x,y) — T^(x,y) is a particular solution of the equation 
(1.1) and T, T p and T^ satisfy the inequalities (2.4). 

By the comment preceding the theorem, 



oc 



T p (x, y) := T(x, y) - T h (x, y) = - £(-!)' f* + j ™ + J ' g) . - (2.23) 

j=o HLo^^ + ^y + ^j 

is well defined and a particular solution of the equation (1.1). 

We also have the inequalities (2.4) by taking the limit on both sides in (2.9) and (2.18). 

Step 5. We prove the uniqueness of T,T p and T^ subject to (2.4). 

Now assume that T ,T h ,T are another mappings satisfying the conclusions in the theorem. 
Since T,T satisfy the equation (1.1), 

T(x + 2mp,y + 2mq) y^ ■ ip(x + jp,y + jq) 

U 2 i=o~ ll P(x + ip,y + iq) j=0 Y]? i=0 V>{x + ip,y + iq) 



and 
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, T'(x + 2mp,y + 2mq) ^ ■ ip(x + jp, y + jq) 



X\i=$ <f{x + ip,y + iq) p£ \{\ =0 ^{x + ip,y + iq) 

for any x,y £ Nk and thus it then follows from (2.4) that 

/ \T(x + 2mp,y + 2mq) -T'(x + 2mp,y + 2mq)\ 

\T(x,y)-T(x,y)\ = 2m-1 , — (2.24) 

lli=o ^{x + ip,y + tq) 



< 



25 AA 1 



En 



Ui=o f(x + ip,y + iq) f^ l= o <P( X + 2m P + iPiV + 2m( l + »?) 
— > as m — > oo. 

This implies the uniqueness of T. Similarly we have the uniqueness of T/j by the same method of 
(2.24). This completes the proof of the theorem. □ 

Remark 2.2 // ip = 0, then T = Th, T p = since e (x, y) = /or any x, y € J%. 

In the next theorem, let two functions ip : N^ — ► (0, oo), ^ : Nfc — ► M satisfy 

oo 3 -. 

,'(x):=ft^±M<0 (2.26) 

j=0 lli=0 <P( X + *P) 

for all x G Nfc. Then we obtain the Hyers-Ulam stability problem for a single variable. 

Theorem 2.3 Suppose that functions f,ip:Nk^R and if satisfy the inequality 

\f{x + p) + ip{x)f{x) + il){x)\ <S 

for all x G Nfc. Then there exist unique functions T,Th,T p : N& — ► K such that T,T p satisfy the 
equation f(x + p) + <p(x)f(x) + ^(x) = 0, 7\ satisfies the equation f(x + p) + (p(x)f(x) = and 
the relations 

\f(x)-T(x)\<8e(x), 
\f( x )-T h (x)\<5e(x) + \e'(x)\, 
\T p (x)\ < \e\x)\, 
T(x) = T h (x)+T p (x) 

hold for all x G N&. 
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3 Applications 

We list some examples of difference equations which are stable by Theorem 2.1 in the sense of Hyers 
and Ulam. 

The function f(x) = f^(lnt) x dt (x G No) is a solution of the nonhomogeneous difference equation 

/(x + l) + (x + l)/(x)-e = 0. (3.1) 

We obtain from Theorem 2.3 the Hyers-Ulam stability for a single variable as follows. 

Corollary 3.1 Suppose that a function f : No — ► M satisfies the inequality 



\f{x + l) + {x + l)f{x)-e\ <5 



(3.2) 



for all x G No- Then there exist unique functions T,Th,T p : No — ► K such that T,T P satisfy the 
equation (3.1), T^ satisfies the equation /(x + 1) + (x + l)/(x) = and the relations 



\f(x) - T(x)\ < (e - 1)5, \f(x) - T h (x)\ < (e - 1)5 + e, 

|T p (x)| < e, T(x) = T h {x) + T p (x) 



(3.3) 



hold for all x G Nq- 



Proof. We apply Theorem 2.3 with <p(x) = x + 1 and V'(x) = — e. For any x G No 

oo j 



yn — x — 



1 



+ 



1 



j=0 i=0 



and 



(-iy(-e) 



^oIll=o x + i + 1 



Thus we lead to the conclusion. 



x + 1 (x+l)(x + 2; 
< e-1 



' 'i-^+ ' 



+ 



(3.4) 



x + 1 V x + 2 (x + 2)(x + 3) 
< e. 



+ • 



(3.5) 



D 



The function /(x) = L t x e l 1 dt is a solution of the nonhomogeneous difference equation 

/(x + l) + (x+l)/(x)-l=0 (xGN ). (3.6) 

We obtain from Theorem 2.3 the Hyers-Ulam stability of (3.6), of which the proof is similar to that 
of Corollary 3.1. 
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Corollary 3.2 Suppose that a function f : Nq — ► M satisfies the inequality 



\f(x + l) + (x + l)f(x)-l\<5 



(3.7) 



/or a// x G No- T/ien i/iere exist unique functions T,T^,Tp : No — ► R suc/i i/tai T, T p satisfy the 
equation (3.6), Th satisfies the equation f{x + 1) + (x + l)/(x) = and i/ie relations 



|/(x) - T(x)| < (e - 1)5, |/(x) - T ft (x)| < (e - 1)5 + 1, 

|T p (x)|<l, T(x) = r h (x) + T p (x) 



(3.8) 



hold for all x G No- 

The function /(x) = L 4+idt (x G N) is a solution of the nonhomogeneous difference equation 



/(x + l) + 5/(x) 



1 



x + 1 



(3.9) 



We obtain from Theorem 2.3 the Hyers-Ulam stability of (3.9) as follows. 

Corollary 3.3 Suppose that a function f : N — > M satisfies the inequality 

1 



/(x+l) + 5/(x) 



x+1 



< 5 



(3.10) 



/or a// x G N. T/ien t/iere exist unique functions T,Th,T p : N — > M snc/i i/iai T, T p satisfy the 
equation (3.9), T/j satisfies the equation f{x + 1) + 5/(x) = and t/ie relations 



\f(x)-T(x)\< 6 -, \f( x )-T h (x)\< 5 - + ± 
\T P (x)\<^, T{x) = T h {x) + T p {x) 



hold for all x G N. 



(3.11) 



Proof. We apply Theorem 2.3 with y(x) = 5 and VK^) = ~+i- For an y ^ G N 

A^l 1 1 1 



j=0 i=Q 



5 5 5^ 



and 



E 

3=0 



-IV 



-1 



rr+j+1 



nL 5 



5 V^ + l 5(x + 2) + '"J - io" 



(3.12) 



(3.13) 
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This leads to the conclusion. 



□ 



The function f(x) = f^° ^ ^ dt (x £ N2) is a solution of the nonhomogeneous difference 
equation 

9 



2r — 3 
/(x+l) + f{x) 



0. 



(3.14) 



Here we denote k := Y^^n =rr — — 

^3-° nj= 2i + l 

stability of (3.14) from Theorem 2.3. 



x x 

for the sake of abbreviation. We obtain the Hyers-Ulam 



Corollary 3.4 Suppose that a function f : N2 



2r — 3 
/(*+!) + f{x 



satisfies the inequality 
< 5 



(3.15) 



for all x € N2. Then there exist unique functions T,Th,T p : N2 — ► K suc/i i/tai T, T p satisfy the 
equation (3.14), Th satisfies the equation fix + 1) -\ — ^F^/(x) = and i/ie relations 



|/(x) - T(x)| < K <5, |/(x) - T fc (x)| < k5 + 9, 

|T p (x)|<9, Tix) = T h ix)+T P ix) 



hold for all x S N2. 

Proof. We apply Theorem 2.3 with ^>(x) 

x + i 



00 j 

En 

j=0 i=0 



"^A and ipix) = —. For any 1GN2 
x(x + 1) 



x 



2x + 2i - 3 



+ 



2x-3 (2x-3)(2x-l) 



+ 



< k 



(3.16) 



(3.17) 



and 



oc 



(-1V-^ 

V ' X+'i 



x+j 



nj 2x+2i-3 
i=0 rr+i 



This yields the conclusion. 



9 



1 



+ 



x(x + 1) 



2x-3V 2x-l (2x-l)(2x + l) 
< 9. 



+ 



(3.18) 



D 



The function /(x) 
equation 



/1 



00 1 



t a: v / 3+6* 



dt (x € N) is a solution of the nonhomogeneous difference 



2r — 1 I 

/(x+l) + /(x)-- = 0. 

X X 



(3.19) 
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Here we denote uj = y~*°° , —j J ' for abbreviation. We obtain from Theorem 2.3 the Hyers-Ulam 

stability as follows. The proof of Corollary 3.5 is similar to that of Corollary 3.4. 

Corollary 3.5 Suppose that a function f : N — ► M. satisfies the inequality 

2r — 1 1 

|/(x + l) + f( x )--\<5 (3.20) 

x x 

for all x G N. Then there exist unique functions T,Th,T p : N — > M suc/i i/iai T, T p satisfy the 
equation (3.19), T^ satisfies the equation f{x + 1) + 2x ~ 1 f(x) = and the relations 

\f(x)-T(x)\<uS, \f(x)-T h (x)\<uS+l, (3.21) 

\T p (x)\<l, T(x) = T h (x) + T p (x) 

hold for all x € N. 
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1. Introduction 

The following two theorems are important in coding theory (see, for instance, 

[1], [4, p. 43-47]): 

Theorem 1.1. (Kraft's Theorem) We have the following: 

(1) If C is an r-ary instantaneous code with codeword lengths li, I2, • • • , l n , then 
these lengths must satisfy Kraft's inequality: 



n 



< 1. 

fc=i 



(2) If the numbers hjh,--- ,ln o,nd r satisfying Kraft's inequality (1.1), then 
there is an instantaneous r-ary code with codeword lengths h,h,- • • J n - 
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Theorem 1.2. (McMillan's Theorem) If C = {ci, C2, • • • , c n } is a uniquely deci- 
pherable r-ary code, then its codeword lengths must satisfy Kraft's inequality (1.1). 

Lei us consider, for an r-ary code C having the codeword lengths l±, 1%, • • • , l n , 
the Kraft numbers 



n 1 

K r (l 1 ,l 2 ,...,l n ) = Y,-r (r>l). 



k=i 



The following estimates of these numbers are given in [1]: 

Theorem 1.3. Let C = {c\,C2,--- ,c n } be an r-ary code having the codeword 
lengths l\, l 2 , • • • J n - Then we have the following inequalities: 

1 n 

— ^[lnH-I^lnr) 2 ] 

k=i 

< K r (li,l 2 , ■ ■ ■ Jn) 

1 ^-^ r li In r + n\nn — nli(lnr) 

~ n lnr ^-— ' r li 

fc=i L 

The equality holds if and only if U = log r n. 

In this paper, we shall give some further improvements of Theorem 1.3 as well 
as some related results. 

II. The Main Results 

First, let us prove the following: 

Theorem 2.1. Let C = (ci,c 2 ,--- ,c n ) be an r-ary code having the codeword 
lengths l\, I2, • • • J n - Then we have the following inequalities: 

nr-™^=i li <K r (l u l 2 ,--- ,l n ) 
2-1) 

The equality holds if and only if h — h — ■ ■ • — In- 
Proof. We shall use the well-known Jensen inequality for convex functions / : J — > 
R, that is, 

(2-2) /gX>)^E/< 



J i 



where / is an interval of R and Xi G / for i = 1, 2, • • • , n. 
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For f(x) = r x , if x% = —1% for all i = 1,2, •• • , n, then we have the first inee- 
quality in (2.1). Note also that the function f(x) = xlog r x is convex since r > 1. 
Thus the inequality (2.2) for Xi = r~ li (i = 1,2, ••• ,n) becomes the following 
inequalities: 



i.e. 



i=l \ i=l / i-\ 



^[Itr-Ai Eil '"" 



rp H 



nfrt J E<=i»- 

and so we have the second inequality in (2.1). This completes the proof. 

As consequences of Theorem 2.1, we have the following corollaries from [1]: 

Corollary 2.2. Let C = (ci,C2,--- , c n ) &e an r-ary code having the codeword 
lengths h,h,- • • J n - If we have the following inequalities: 



1 n 
(2.3) -J><log r 



then C is not uniquely decipherable. 

Corollary 2.3. If the real numbers r and U (i = 1, 2, • • • , n) satisfy the following 
inequality: 

(2-4) ^= 1 * - >log r n, 

t/ien there is an instantaneous r-ary code with codeword lengths l±, I2, • • • ,l n - 

An extension of Theorem 1.3 is given by the following: 

Theorem 2.4. Let C = (ci,C2,--- , c n ) be an r-ary code having the codeword 
lengths l\, I2, • • • J n - For arbitrary y,z 6R, we have the following inequalities: 

(2.5) F(y)<K r (hJ 2 ,--- ,l n )<G(z), 
where 

(2.6) F(y) = nr y — r y lnr 2_, h + n V 



(2.7) G(z) = nr z -J2 r ~ h (h + z)lnr. 
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n 

Ci 



The strongest inequalities of this form (2.5), that is, 

F(y) = maxF(y) < K r (l u l 2 ,--- ,l n ) 

(2.8) y 
V ' <G{z) = minG(z), 

z 

are obtained for 

(2.9) y = --^2l u ~z = log r (-J2r 

i=l \ i=l 

Proof. Let us start, as in N. M. Dragomir et al. [1], from the following inequalities 
for a differential) le and convex function / : D(f) C R — *■ M, 

f'(y)(x -y)< fix) - f(y) < f'(x)(x - y) 

for all x,y G D(f) (: the domain of /). Then, letting f(x) = r x in the first 
inequality and f(x) = r x , y = z, in the second inequality, then we have 

(2.10) r y +r v (x — y) lnr < r x < r 2 +r x (x — z) lnr. 

Also setting x = —k (i — 1, 2, • • • , n) and adding all such inequalities, we have 

n 

(2.11) F{y)<Y j r- li <G{z). 

i=l 

Further we have also 



F'{y) = -r y (lnr) 2 \S^k + ny 
F"(y) — (lnr) < — r y lnr 2_, h + ny \ — r 



i, \ 

^ y 



.i=i 



From F'(y) = 0, it follows that y = y, while F"(y) = — r y (lnr) 2 < and so we 
have the maximum. Similarly, we have also 



n 

h 



G' (y) = nr z In r — In r \_, r 

i=l 

G"(y) = nr z (lnr) 2 . 
From G'{z) = 0, it follows that the minimum value is obtained for 

n 

nr z = y r~ % 

i=l 

that is, for z = z. This completes the proof. 
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Remark 2.1. Note that the first inequality in (2.8) is equivalent to the first 
inequality in (2.1) and so the following improvement of the first inequality in (2.5) 
is valid: 



(2.12) 



F(y)<nr-*Z? = iU < K r (h,l 2 ,- ■ ■ ,l n ). 



Moreover, this result is a simple consequence of the first inequality in (2.1) and 
the first inequality in (2.10). Indeed, setting 



x 



1 n 



i=l 



in the first inequality in (2.10), we have 



r v _ r v 



I — 2_. /i + y J In r < r n ^ ■ 



1 /• 



which is equivalent to the first inequality in (2.12). 

For y = z = log r ( ^ ), (2.12) together with the second inequality in (2.5) gives 
the following inequalities: 



(2.13) 



n In r ^— ' L 


'<" ; f-iV 111 ' ; J 




_ A. y* n i . 




< K r (li,l 2 , ■■ ■ Jn) 




1 U 

< T 

nlnr ^— ' 


r Zi lnr + nlnn - 


- n/i(lnr) 2 


^.tg 


_ 



Remark 2.2. The second inequality in (2.8), after some elementary transforma- 
tions, gives the second inequality in (2.1). The upper bounds for K r (li, 1%, ■ ■ ■ ,l n ) 
are nontrivial if they are < 1. For example, for (2.1), we should have 

(2.14) m .-(£?=i'*r-'<)/(£?=ir- , *)< lj 

while, for (2.13), we should have 



(2.15) 



n 

— Y 



nlnr 



i=l 



r li In r + n In n — n/^ (In r) 



< 1. 



The first condition, (2.14) and (2.15) give (2.4), which is the conclusion of Corollary 
2.3. 

Since the expression on the right-hand side of (2.4) is not simple one, we shall 
give some simple results. First, we prove the following: 
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Theorem 2.5. Let the assumptions of Theorem 2.1 be satisfied and let 

(2.16) m<k<M (i= 1,2,- •• ,n). 

Then we have the following inequality: 

(2.17) K r {h,l 2 ,--- ,Z n )<nr"^^=i\ 
where 

Vme- m + \lMe~ M 



(2.18) K 



Vme- M + VMe" 



Proof. Note that the following result is valid (see, for instance, [2, p. 307]): 
If < a < ai < A and < b < bi < B for all i = 1, 2, • • • , n, then we have 

(2-19) J- < " E ^f* =- < K 2 , 

where 



^ = ^M (>d. 



Setting a,i = li,bi = r li {% — 1, 2, • • • , n), a = m, ^4 = M, b = r M and B = r~ 
then the first inequality in (2.19) gives the following inequality: 

V n 1r~ l i 1 1 J^ 



Ei.r-'' "if 2 



■/) 



1=1 



Therefore, from the second inequality in (2.1), we obtain (2.17). This completes 
the proof. 

Corollary 2.6. If the real numbers r and li (i = 1, 2, • • • ,n) such that < m < 
li < M and r > 1 satisfy the following inequality: 



1 n 
(2.21) -V/,>K 2 log 



i=l 



where K is defined by (2.18), then there is an instantaneous r-ary code with code- 
word lengths h,h,- • ■ , In- 
Proof. The inequalities (2.20) and (2.21) gives the conclusion (2.4). 
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Theorem 2.7. Let the assumptions of Theorem 2.5 be satisfied. Then we have 
the following inequalities: 

(2.22) K r (/i,/ 2 ,--- »«n) < ^ r- H — r 

M — m M — m 

Proof. The following result for a convex function / : J — > R is valid [3, p. 98]: If 
m < Xi < M for m, M E I and i = 1, 2, • • • , n, then we have 

(2-23) I £ /(*,) < M "J E " =1 g V (m) + " E y gi " m /(M). 

n ^^ M — m M — m 

i=l 

Letting f(x) = r x ,Xi = —li (i = 1,2, ••• ,n), M — ► —to and to — ► — M, then we 
obtain (2.22). This completes the proof. 

Corollary 2.8. // t/ie rea/ numbers r and li (i — 1, 2, • • • , n) siic/i t/iat < to < 
U < M and r > 1 satisfy the following inequality: 

lA Mr" m - TOr" m _ 1 M - to 

V * ' „, / / l j, — m y — M rr. y — m y — M ' 

1=1 

t/ien there is an instantaneous r-ary code with codeword lengths h,h,- • • ,l n - 
Proof. It is easy to see that (2.24) is equivalent to the following: 

nM-YT l =xU r - m + Y,?=ik-nm r _ M < 1 
M — to M — m ~ 

and so (2.22) gives 

K r (li,l 2 , ■ ■ ■ ,l n ) < 1. 

Applying Kraft's theorem, we get the above result. This completes the proof. 

Theorem 2.9. Let the assumptions of Theorem 2.1 be satisfied. Then we have 
the following inequality: 



Kr(h, h, • • • Jr 



n n 



<nr nT,7=i l i - l n rS^Ur h -\ V" U V" r h 

i=l i=l i=l 

If m < U < M for all % = 1, 2, • • • , n, t/ien twe /lave £/ie following inequality: 



K r (h, h, ■ • ■ ,l n ) 

(2.26) _j_v"™ ; l nr 
<^ <fifr n ^^—i % -4- 

~~ ' n 



r 2 



(M-m)(r" m -r- M )- 
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Proof. The inequality (2.25) is a special case of the second inequality in (2.5) for 
z = — ^ Y^H-i h- For ^ ne P ro °f °f (2-26), we need the following results ([2, p. 299], 
[3, p. 206]): Let a and b be n-tuples such that u < cii < U and v < bi < V for all 
i — 1, 2, • • • , n and real numbers n, v, U, V. Then we have the following inequality: 

,2- 



(2.27) 



IV j / V IV 

J2 a * b *--J2 a *J2 b * 



< 



1 



■/? 



n 



(U-u)(V-v). 



Note that, for a^ = U and bi = r Ci (i — 1, 2, • • • , n), we have £/ = M, u = m 
Y = r — m r = r — M 
{/i} and {i; *}, we have 



r = r ■" and, by Cebysev's inequality, for oppositely ordered n-tuples 



-L- 



Therefore, (2.27) becomes the following inequality: 



n n 



i£<.£>-<.-5> 



i=l 



n n 



j / ^ lb IV 

^E'-E^-E'* 



< 



» 



n 



(M -m)(r 



-h 



-m _ r -M^ 



and, from (2.25), we have (2.26). This completes the proof. 

Corollary 2.10. If the real numbers r > 1 and U (i = 1,2, 
following inequality: 

n 



, n) satisfy the 



lnr 



n n 



nr 



E"=i^ < i + lnr V" lnr li X^*X^ 



or i/ i/ie rea/ numbers k (i — 1, 2, • • • , n) are siic/i £/ia£ m < k < M and satisfy the 
following inequality: 

lnr 



nr 






r 2n 



n 



(M -m)(r 



-m _ -M\ 



then there is an instantaneous r-ary code with codeword lengths l±, I2, • ■ • , l n - 
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1 Introduction 

In this paper we consider approximations to the solution of a system of ordinary differential equations 
(ODE) in presence of possibly infinitely many state dependent impulses on the righthand side [7], 

[8]: 

f(x(t),t)+ £ a {x{Tr))S{t-T ), te [0,71, 

3 eJ (1) 




where T > is a real number , x° G E is the initial condition, where E = M. d with d a positive 
integer, x : [0, T] — > E is a vector function , / : E x [0, T] — ► E is a given map, 8 (.) is the Dirac 
delta distribution at 0, J C N* = {1, 2, 3, ...} is a countable set of indices, oij : E — ► E is a given 
map for all j € J and {tj}j<=j is a strictly increasing sequence on ]0,T]. 

In [6], the variational formulation is applied to obtain existence and uniqueness of the solution of 
(1). Then, using a Galerkin method with piecewise polynomial functions of degree K , it is proved 
that the rate of convergence in the L 2 -norm is h 1 / 2 and the rate of convergence at the nodes is h 1 . 
In the state independent impulses case the rate of convergence at the nodes grows up to h K+2 [3] . 
Moreover, when the solution is in H K+1 (0, T; E), the rate of convergence in the L 2 -norm and at the 
nodes become respectively h K+1 and h 2K+2 [2] , [5]. 

In this paper we use a simple transformation that allows us to show that the problem is equivalent 
to an ordinary differential equation without impulse with its solution in ^(O^-^E). This trans- 
formation is numerically important. Indeed, to approximate the solution of (1), we just have to 
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... approximate the solution of the equivalent problem which is regular. We then obtain better rates of 
convergence. Moreover we obtain an approximation of the solution which looks like the solution of 
(1). Both, the solution and its approximation, are in BV(0, T; E), the space of functions of bounded 
variation on [0, T] . 

In section 2, we present the equivalent problem without impulses (1). Its variational formulation is 
given in section 3, and we show existence and uniqueness of the solution under weaker assumptions 
than those required in [6]. In section 4, we consider the approximation problem by using piecewise 
polynomial functions. Then, in section 5, we show that the L 2 and the nodal convergence rates are 
h . Moreover, in the state-independent impulse case, we show that the nodal convergence rate is 
h 2 for nonlinear systems and h K+2 for linear systems. Numerical implementation is considered in 
section 6 and examples are given in section 7. 



2 Transformation of the problem into an ordinary differential 
equation without impulse 

A solution x of (1) is a function of bounded variation which is the sum of two parts: x(t) = 
x n(t) + Xj(t), a continuous part Xr in H l (Q,T\E) (the Sobolev space of L 2 -function with L 2 - 
derivative) Xr{€) = x° + L /(x(t), r)dr, and a jump part xj in L°°(0, T; E) (the space of essentially 
bounded functions) where Xj(t) = "}2 aj(x(T~))x[ T -,+cc)(t)- 

We introduce two families of functions: {0j}j^j and {(3j}je.j from C(0, T; E) to E defined by 

( 6 1 (y)=y(T 1 ), 

< e j (y) = y(T j ) + Ea i °0i(y) forj>l, 

I i<3 

and (3j = ctj o 6j for j e J. If x is the solution of ( 1 ) , then a j (x(t~ ) ) = (3j(x), and xr is the solution 
of the following ODE 

x R (t)=F{x R ,t), t£[0,T], 

x R (0) - x°, (Z) 

where F : C(0,T;E) x [0,T] — » E and F(y,t) = f(y(t) + £ P J (y) X [r 1 ,+oo) (*),«)• 

3 Weak Formulation 

Using the Galcrkin method we approach the solution of (2) by a piecewise polynomial function u^ . 

The values Uhijj) approach xr(tj) (j € J). Then we show that Uh + Yl Pj( u h)X[T- +oo) converges 

jeJ 
to the solution of (1). 

The (global) weak variational formulation of (2) is 

find (x R ,X) e ff^O.T;^) x E such that 

X ■ v(T) - / T xr(t) ■ v(r)dT - x° ■ V (0) + tf F{x R , r) • v(r)dT, (3) 

for all v e ^(O.T;^). 

Theorem 3.1. Problems (2) and (3) are equivalent. 

Proof : Sec [3, Theorem 3.3]. □ 

A partition of the interval [0, T] of size h is characterized by an integer N > 1 and a sequence 
{tn}n=o °f re& l numbers such that = to < ■ ■ ■ < t n < • • • < £jv = T, with h = max{£„ — t n -\ : n = 
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1, 2, ... , N}. Subintervals will be denoted by /„ =] t n -i,t n [ , n = 1, . . . , N. A family of partitions ... 
indexed by h is said to be regular if there exists a constant c > such that as h goes to zero 
ch < t n — t n -i < h, n = 1, . . . , N. Throughout this paper we shall only consider regular families of 
partitions. 

n 

For any n = 1,...,N, let U n = E n+1 x Jl H 1 (Ik;E) and for u n = (Uq, ..., U n ; m, ...,«„) G U n we 

fe=i 
define 

. r .x_ f 17* ift = * fc (fc = 0,...,n), 

MnW 'l «*(*) ift€]*fc_i,t fc [ (fc = l,...,n). 

To obtain a piecewise polynomial approximation of the solution, we consider the following variational 
formulation 

' If U = x°, and (u k , U k ) G H l {I k \ E) x E is given 
for fc = 1,2, ....n — 1, find (« tt , [/„) G H x {I n -,E) x E such that 
U n -v n (t n )-fj u„ (t) • v n (r)dT = U n -\ ■ v n (t n -i) + Jj F(u n ,r) ■ v n (r)dT, 

{ for all v n eH?(I n ;E). 

If this problem has a solution, we obtain, as in the proof of Theorem 3.1, U n = u n (t n ) and U n -i = 
u„(i n _i). Thus, we can easily show that u n {.) is continuous and belongs to iJ 1 (0,<„;_E). 

Theorem 3.2. Problems (3) and (4) are equivalent. 

Proof : Sec [3,Thcorem 3.3]. □ 

Let us define 

J n = {j G J : Tj G ]*„_i, t„]} , n = 1, ..., N, 

j£={j€J:T j €]0,t n }}, n = l,...,N 

and 

lco(E J ) - {? G £ J :|| C ||oo= *«p{|&| : J e J} < +^}- 

By the same way we define lao(E Jn ) and loo(E J ™). 

To show existence and uniqueness of the solution to (4) we use the next lemmas [8] . To simplify 

+oo 

we suppose J = N*, however our results remain true for any subset J of N*. Let A = J^ Aj, and set 

n n 

Q 7j = 1 if n < m and ^2 ji = if n < m. 

i—'m i—'tn 

Lemma 3.1. Let {Xj}j£j be a sequence of positive real numbers , then 

+oo i— 1 



l + ^A,[](l + A fe )< e A . D 



j=i fe=i 

Let us remark that the functions 0j(.) are well defined only if the values y(rj),j G J, are defined. 
For a sequence £ = {£j}jej we can define 

i ^(0 = ^ + Eai°^(0 forj>l, 

with, as previously defined, /?j = a-,- o 0j for j G J. 
Lemma 3.2. For any £ x and £ 2 in looiE-'), we have 

\m l )-m 2 )\ < A j e A max{|^-^|:fc = l,...,i} 

< A J e A ||e 1 -e 2 ||co. □ 
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Remark 3.1. If x\,x 2 £ C(0,T;E) and £* = x k (Tj) (k = 1,2), i/ien 

IA(^)-^(^)I < A J e A max{|xi(r fe )-o ; 2(r fe )|:fc = l,...,j} 

D 





< Aj-e A | xi - x 2 oo- 




Lemma 3.3. 


We /mue J2 l&(0)| < e A E |«j(0)|. 


a 


Lemma 3.4. 


Let I = [a, 13]. 




(%) TVie map 








H\I-E) — » L 2 (/;£) 


X £ 




u i — ► {v,v(0)) 





is an isomorphism. 

(ii) Let b be an arbitrary element of H (I,E)*, the variational problem 

find (u, U) £ L 2 (I; E) x E such that 

U.v{(3) - f a u(t) ■ v{t)dt = b(v) (5) 

for all v£ H^E) 

has a unique solution. 

Proof : Sec [10]. □ 

Theorem 3.3. Let x° £ E, 

(a) assume the map f : E x [0, T] — ► E verifies the following conditions : 
(i) for all x £ E, the map t — ► f(x,t) is ( Lebesgue ) measurable, 

(ii) there exists q £ L 2 (0,T;R) such that, for any X\,x 2 £ L (0,T;E) 

\f{x 1 (t),t)-f(x 2 {t),t) \<q{t) \ Xl (t) - x 2 (t) \ a.e., 

(Hi) the map t ► /(0, t) belongs to L 2 (0,T;E); 

(b) suppose also the functions ctj (j £ J) verify the following conditions : 

M E \<*j(o)\ < +oo, 

jeJ 

(v) there exists a sequence of positive real numbers {Xj}jej such that ^ \j < +oo and for all j £ J 

and xi,x 2 £ E 

| Ctj(xi) - Oij(x 2 ) \< Xj | X\ — x 2 I • 

Then 

(1) the solution {U\, ..., Ujf, U\, ..., Un) £ Un to (4) exists and is unique, moreover, xr — Ujv is the 
unique solution in H 1 (0,T;E) to the system (2), and 

(2) the solution to (1) exists and is unique in BV(0,T; E), moreover, we have 
x(t) = u N (t) + J2 ft( u w)X[r j ,+oo)(0- 

Proof : Similar to the proof of Theorem 3.5 in [6]. See also [11]. □ 

Let us consider the Hilbert spaces U and V defined by 

N N 

U = E N+1 x YlH^I^E) = U N and V = E x J[ H 1 (J„ ; E) =V N . 

n—l n—1 

In order to apply a Galcrkin method we consider the following equivalent problem which is a global 
mesh-dependent variational formulation 

find u = (Uq, U\, ..., Un] Ui, •••, Ujv) £ U such that 

N N 

Uq(V - Vi(t )) + J2 U n (v n (t n ) - v n+1 (t n )) + U N ■ v N (t N ) - J2 Ii u„(t) ■ v n (r)dT 



y 



n=l ' ( 6 ) 



X°-V + J2lJ In F(u n ,T)-V n (T)dT}, 



n=l 



for all v = (Vq, Vi, ..., v^) £ V. 
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4 Galerkin Approximation 

Let us define the subspaces U and V of W and V by 

■iih [ u h = {U ,Ui,..,U N ;ui,...,u N )<EU 1 ,. 

1 such that u n e T K (I n ;E) for n = 1, ...TV J ' l J 



and 

such that v„ G V K+1 (I n ;E); for n= 1,...,7V 
where T jK (I n \ E) is the set of polynomials of degree at most K on I n . For each n = 1, ..., TV we set 



yh = j V h = (V ,V 1 ,V 2 ,...,V N ) &V i ,, 



M^B^xH^l^), 



fe=l 

h 



and for each u„ = {U$, ..., [/„; ui, ..., u n ) G W„ we define 

J [/ fe if * = *fc (fc = 0,...,n), - , 

MtJ_ i#) ift€]*fc_i,t fc [ (* = l,..,n). W 

The discrete problem associated to (6) is 

find u h = (Uo, ..., [/jv; Wi, ..., Ujv) G W ft such that Uq — x° and for n = 1, ...,7V 
£4 • v„(t„) - Jj u„(r) • v n (T)d,T = f7„_i • w„(t„-i) + J^ F(u n ,r) ■ u„(r)rfr, (10) 

for all u„ eV^ +1 {I n ;E). 

We will show that (10) has a unique solution for /i small enough. Similar to Lemma 3.4 and Theorem 
3.3 we have 

Lemma 4.1. Let I n = [t n -i,t n ]. 
(i) The map 

r K+1 (I n ;E) — V K (I n ;E)xE 

V I ► (-V,v(t n )) 

is an isomorphism. 

(ii) Let b an arbitrary element ofT 1 (I n ; E)* , the variational problem 

find (u, U) G V K (I n ; E) x E such that 

U.v(t n ) - J u(t).v(t)dt = b(v) for all v G V K+1 (I n ;E), 

has a unique solution. □ 

Theorem 4.1. Let the assumptions of Theorem 3.3 be verified, then the variational problem (10) 
has a unique solution. 

Proof : Let Uo = x° and, on each interval /„, we suppose that {v,k, Uk), k = 1, ...,n — 1, are given. 
We look for (u„, U n ) G T K (L n ;E) x £ such that 

U n ■ v n (t n ) - Jj u n (r) ■ v n (T)dT = U n -! ■ v n (t n -i) + Jj F(u„,t) ■ v n {r)dr 
for v n eV K+1 \l n ;E). 

For all u„ G V K (I n ;E), the map u — ► 6„(w;u„) = U n -i.v(t n -i) + f In F(u n ,T) ■ v(r)dr is 
well defined, linear and continuous on V K (I n ',E), thus b n (.\u n ) G T jK+l (I n \E)* . Let (u°,£/°) G 
V K (I n ; E) x E be arbitrary and fixed. We construct the sequence {(u™, U™)}^^ as follows : given 
(u™, U™) G -P K (/„; E) x £, (m™ +1 , t/™ +1 ) G V K (I n ; E) x E is the solution of the problem 

U™ +1 .v n (t n ) - / <* +1 (t) • ti n (r)dr = &„(«„;<), ( n ) 
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for all v n G T K+1 (I„;E). We show that {w™}+~ is a Cauchy sequence in V K (I n ;E). We have 



|f(*c+\t)-f«»| 

< q(r) | u™ +1 (t) - < l (r) | + | E (/W +1 )X[t,, + oo)M - ^(^)x [r ,, +00 )(r)) | 

For w n satisfying u n (t) = — (w™ +1 — w™) (£) in 7 n , v„(£ n ) = 0, and using Lemma 4.1 and (11) , we 
show that 

nc +1 - elk™ < /»y 2 n«iio,n [nc - c-'iio,™ + ^ /2 Ae A nr - r-'iu.n 

with C° = u n ( Tj ),j e J„, and ||C" " C^L^ - «*P{Kfa) - u^fa)!; j e ■£}• For all 
j £ ^rT-i we have w^(tj) = w™ _1 (tj) since Uk : k — l,...,n— 1, are fixed. By using the following 
inequalities [2] 

II/-771 /-m— 1\\ <-- l — 1/2 rll /~m /-m— 1\\ _i_ ^ ||/° m /-m— IN 1 

II' ' lloo.n — n Ml' ' 110, n II' ' lll,n J 

>m /-m— 111 ^- -L — 1 II/-TI /-m— IN 



C m -C L <c/i" i C m -C 

I ^ =• II 1,71 — n II ^ ^ 



10, n 



we obtain ||m™ +1 - u™\\ n < chl/ 2 \\q\\ 0n \\u™ ~ u™ 1 \\ 0n - If h n is small enough then {w™}+= is 
a Cauchy sequence in V K (I n ;E). Then {u™}^^ converge to u„ G ? K (/„; £7). By continuity we 
obtain — J". u n (r) ■ v n (r)dT = b n (v n ;u n ) for all v n G H (I n ;E), such that v n (t n ) = 0. If in (11) 
we take v n (i) =Ve£,wc obtain C/™ +1 .V = 6„(V; <*) for all V in E. Then {i7™}+~ converge to 
U n £ E and [/„.¥ = b n (V;u n ) for all Ve£ Let w n (t) = v n (t) — v n (t n ) with v n G i? 1 (/ n ;S), wc 
have — J, u m (t) • w n {r)dT = b n (w n ; u n ), since w n {t n ) = 0, and we obtain 

w„(t) • v n (r)dT = b n (v n ;u n ) - b n (v n (t n ); u n ) = &„(«„; u„) - U n .v n (t n ). 

Then (u„, [/„) is a solution of (10) on /„. D 

5 L 2 and nodal errors 

5.1 Convergence results 
Theorem 5.1. Let x G H K+1 (0,T;E), then 

Inf{\\x - w„|| ,„ : u n G T K (I n ; E)} < ch k+1 \\x K+1 \\o, n , (12) 

Inf{\\x-u h \\ :u h \ In =u n GP K (I n ;E)}<ch k+1 \\x K+ %. □ (13) 

Let us set e 0) „ = ||«„ - »_r||o,«, eo = ||«n - ErIIoj e 0i „ = ||u„ - jcrHo,™, E n = \U n - x R (t n )\, 
ei,„ = ||w„ - a;/j||i,n) ei,„ = ||w„ - a;je||i,„, e^^ = sup{|m„(t,) - x fl (Tj)| : j G J„} and 
eoo,n< = swp{|w„(rj) - £r(t,)| : j G J^} where u n G V K {I n ;E). 

Theorem 5.2. Suppose the assumptions of Theorem 3.3 are verified. If h is small enough we have 



\\u N -xn(t)\\ <ch x 
max{\U n -x R (t n )\ ;n = 0,1,.., AT} < ch l 



,.(i) 



(14) 
(15) 
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Proof : We use a proof by induction. Suppose that Ek = |£4 — Xit(tk)\ = 0(h), k = 1, ...,n — 1, 
and || tt — £ij||o,[o,i n _i] = 0(h). From the following problems on /„ 

find x n £ H l (I n ;E), with xi, ...,x„_i already known x n — xn\i n , such that 
XR{t n ) ■ v n (t n ) - J Jn x r {t) ■ v n (i~)dT = a;ij(i„_i) • i>„(i„-i) + J Jn F(x r ,t) ■ v n (r)dT, 
for all w„ £ ff 1 (!„;£) 



and 



we obtain 



find «„ s V (I n ', E) with wi, ..., u„_i already known, such that 
U n -v n (t n )-fj u„(t) ■ v n (T)dr = [/■„_! • «„(<„_ 1) + Jj F(u n ,r) ■ v n (r)dr, 
ioiSi\\v n £V^ +1 (I n ;E) 7 



[U„ - X R (t n )} ■ V n (t n ) - J In (Un(T) - Xr(t)) ■ V n (r)dT = [J7 n _i - £r(£„_i)] ■ U„(t re _l) 

+ ]/ [-^(Mn, r) - F(a;i{, r)] • v n (r)dT. 



(16) 



Take w n = U n — xn(t n ) in this equation. Since e^^ < h n [erj, n + ^n 6 !.™], see [2], we obtain the 
following estimates 

E„ < E n _! + (1 + A n e An )\\q\\ , n eo,n + h. n e An h n ei, n \\q\\o tn + & n e A ' n \\q\\o,nh 1/2 e 00tn _ 1 <. 

Moreover, ei, n < ||w n - w„||i,„ + ei,„ and ||u„ -u„||i,„ < cft" 1 ^ - w„||o, n - Then we obtain 

E n < -E„-i + (l+A n e A ")||(7||o,„eo ;Il +A n e A ' l /i^{c^ 1 ||?i„-u„||o,„+ei ! „ + /i _1/2 e 00:n _ 1 <}||g||o,„. (17) 
By taking v n such that v„(t) = — (u„ — u n )(r), t £ I n , and v n (t n ) = in (16), we obtain 

||«n-«n||o,n < J^Z" {/^ /2 -E„-l + (1 + 0„)e O ,„ + ^ /2 |k||o,nA„e A "e li „} 

1 



1-. 



-/in|kllo,nA„e "e oon _ 1 <, 



(18) 



1/2 



where #„ = h n ||g||o,n(l + (1 + c )A„e "). By using this inequality for the L -error, we obtain 



eo.n < 



1 



1- 



hl/ 2 E n _ 1 + 2e ,„ + h n /2 \\q\\ , n A n e A "ei tn + /i^||g|| ,„A„e A "e oo i; 



By the same way for the nodal estimates, we use (18), then (17) becomes 
l + ||g||o,nA T ,e A »ft 1 /2 



E n < E n . 



+eo, n 



-ei„ 



(1 + A„e A ")||g||o.n + T^TrdMMl + A„e A ") + || g || ,„A n e A ") 

1 — f n 



/in||g||o,nA„e A " + ^ /2 ||g|| ,„A„e A " 



lkl|o,n(l + A„e A ») + c|| g || , rl A„e A " 



1-0. 



l + jf^a + AneA-J^HgHo, 



+eoo,n-i<A„e A ™ft,y 2 ||g|| ,, 
For h small enough, the coefficient of e^ „_i< can be rearranged as follows 



K n e^h}J 2 \\q\\ , 



1 + rVi 1 + A„e A ")/ l y 2 |M| ,„ 

J- — f n 



< A„ e A ^ r 1 / 2 || g || .„(l + -^— 0„), 

1 — Cn 



< 



c'^ 



Approximation of State-Dependent... 



424 



The coefficient of eo,„ is 

(1 + A„^) Mk „ | (1 _ K) + (1 + 8n) + c(1 + eu) _^l_} < (1 + A .^. )| , |0|| (3 + ft.) 

and the coefficient of ei in is 

/i„||g||o, n A„e A 



l + A„e A - 



1 



1 - 0„ + ftX /2 + c/iy 2 ||g|| ,„A„ e / 



< 



fenllgHo.nAne 7 ' 
1-ft,. 



l + c/i I 1 / 2 ||q|| ,„A„e / 



Ifwetakc a „ - i±|a, /?„ - IMIo,n[i+Ane A "c] (3+ ^ 7n = _J_ r^y^ „A ne A "(l+|| <z || ,„^ /2 A„ e A " C )} 
and <$ n = -^J 1 , we obtain E n < a n E n ^i + /3 n eo, n + 7 n ei. n + ^n^oo n _i< • Step by step we show that 

n n n 

E n < KE + K J2 Pn-jeo.n-j + K E Jn- 3 ei.n-j + K E ^n-j&oo n-j<i where K is defined by 

3=0 3=0 





3=0 


K = 


n 2 /,— 

-- n ^ < e * =1 



»=1 



_ i+e„ _ 



26>„ 



If ft. is small enough so that ft < 1/2 for all I in 



4E»* 



the set {1, ..., N}, then if < e i=1 < +oo. We also have E 0n-j&o,n-j < (E Pn-iY^^o where 



3=0 



3=0 



eo 



= (E eo.ra-j) 1 ^ 2 - By the same way, if we take ei = (E e^n-j) 1 ^ we obtain E 7„_jei )n _j < 



3=0 



3=0 



3=0 



(E 7n-,) 1/2 ei- We show that 



3=0 



(3 + ft) 



E^-j- = E^^ max {7i3^ :0 ^^ ri }EiK™( 1+A ™ eAmc ) 2 - 



i=o 



m=0 



m=0 



If /lis small enough, such that ft < 1/2, we have E Pn-j — c(l+Aec) 2 ||g||§. Since l+/i„ A„e "c||g||o,n < 



1/2 A A 



j=o 



1 + 6>„, then 7„ < h n T ^||q|| ,„A„e A "(l + 6»„), and we obtai 



ail) 



E^^iMiUAAfft 2 . 

3=0 

71, TL Ti Ti 

We also have E <*n-j = E 4 = E f^ < 2c ' E ftfc < 2/ic / ||g|| (l + cAe A ). Finally 
j=o fe=o fe=o fc fc=0 

#n < #£o + Kc(l + Ae A c)|| g ||oe + c|| g || (Ae A c)/i 1 e 1 + K2c'\\q\\ Q h 1 {l + cAe A )e 00 ^„_ J -< . 
We have e^^ < ft™ [e ,„ + h n e^ n ], moreover e 1 , n < ei, n + cft- 1 ^ - w„||o, n , then 
eoo,« < /i~ 1/2 e ,„ + /iy 2 ei,„ + c/i~ 1/2 ||u„ -u n \\o,n, 

< 2h- 1 / 2 e , n + /j-VaiLM |/ l V2^_ 1 + (1 + ^^ + / l 3/ 2 || g || , l A„e A "e 1 ,„ 

+ ft-r l ||<7||o,nA„e A "e oon _ 1 < L 



< 



1 — tf „. 



On the other hand 



eo.n < 



1- 



/i 1/2 £„_i + 2e ,„ + ^ /2 ||g||o,nA„e A "ei ;Il + h n \\q\\ . n K n e Kn e^^^ 
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e < C 



max 



E n + 2e + h 3 ' 2 max{||g|| ,„A„e A " T n = 1, ..., N} ei 



+/i ? 1 / 2 niax{||g|| 0; „A„e A ",n= 1, ..., Nje^ n _ x < 

Let us go back to our proof by induction. For n — 1 we have Eq = Uq — x° = 0, then E\ = 0(h), 
therefore the result is true for E\ and erj,i. Let us suppose that it is true up to n — 1 and we 
show that it is true for n. Since E\ = 0(h) and e^i = 0(h) then e^i < 0(h 1/>2 ). Consequently 
eoo,2 < 0(/i 1/2 ). Finally e^i = 0(h x l 2 ) for all I < n and 

E n < C(E + O(h 1 ) + O(h 3 / 2 )), 

e <[Eo + 0(h n ) + 0(h^ 2 )}. D 
Theorem 5.3. Under the assumptions of Theorem 3.3 and if h is small enough, we have 

\\un + ^2 I3j(un)X[t 3 , +oo) ~x\\o < ch 1 llx^Ho, 



j£J 



,.(D 



max{|C/„ + y^/?j(Miv)X[ r ,., +00 )ftn) -x(£„)|;ra = 0,1, ..,N} < ch 1 

where x is the solution of (1). 

Proof : We have x(t) = x R (t) + £ 0j(x R )X[r :j ,+oo)(tn), then 

jeJ 

\\ U N + ^2Pj(un)X[t j ,+oo) ~x\\ < \\U N -Xr\\ q + \\^2[(3j(UN) ~/3j(a;fl)]X[r ;( ,+cx))l|o- 

jeJ jeJ 

Using the definition of the norm ||.||o, we obtain 

( N 2 1 

|| ^2[/3j(un) - Pj(xr)}X[tj,+°o)\\ = \^2\\Y1 ^o^n) - /3 J (a;i?)]X[r 3 ,+oo)||o i „ \ 



(19) 
(20) 



1/2 



jeJ n 

N 

< e A eoo,w<<! Z!ll Z! A illL 



1/2 



n=l je./„ 
< ^ 1 / 2 Ae A e COiW <. 

From the results of the preceding theorem, we obtain (19). Finally (20) follows. 



□ 



5.2 Superconvergence results for nonlinear systems 

In this paragraph we suppose the impulses are state-independent, ctj(x(r~)) = ctj,j G J. In this 
case [3j(xr) — ctj,j £ J, and the assumption (b) of Theorem 3.3 becomes ^ \otj\ < +oo. With this 

condition, (10) reduces to the following problem 

find u = (Uo, ..., Un; u\, ..., un) G U h such that Uq = x° and for n — 1, ...,N 
U„ ■ v n (t n ) - J In u n (r) ■ v n (r)dT 

= Un-\ ■ V n (t n -l) + fj f(Un(r) + E (XjXfa.+oo)^)^) ■ v n (r)dT, ( 21 ) 



3£J 



for all v n €V K+1 (I n ;E) 



For nonlinear systems, if we use the results of Theorem 5.2 and consider additional conditions on 
F, we obtain higher rates of convergence. 
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Theorem 5.4. If A = (aij)i<ij<d, where a^ G L 2 (0,T;R) and f G L 2 (0,T;E) then (i) the 
unique solution w of the system, w{t) — A(t)w(t) + f(t), w(0) = a, is in H (0,T; E), (ii) the map 
(x , /) i— ► w from E x L 2 {0,T;E) to iJ 1 (0,T;i?) is linear and continuous, and (in) there exists a 
constant c such that 

Nli<c{H + ||/||o}. 



Proof : See [4]. 



□ 



Lemma 5.1. Let x be a function defined from [0, T] to E, and let x^ be its Lagrange interpolation 
of degree on each subinterval I n . If x € H (I n ;E), we have 

\\Xn - x\\o,n < C^H^Ho.n < ch 1 \\x\\ lt „ 

and 

Wx 1 - x\\ Q < ch'Wx'Wo < ch'Wx^. 



Proof: See [1]. 



D 



Theorem 5.5. Suppose the solution of (2) is in -ff 1 (0, T; E), the solution of (6) exists and is unique, 
and the following conditions hold: 

(i) the matrix A(t) = F x (x, t), such that a,ij(i) = ^r(x, t), exists and its columns are in L 2 (0, T; E), 
(ii) there exists a neighborhood V of the origin of E and a constant B > such that for any t G [0, T] 
and y € xr{€) + V 

\F{y,t)~F{x R ,t)~A{t){y-x R {t))\<B\y-x R {t)\ 2 . 

There exist a constant c independent of h such that for h small enough 

max{|C/„ -x R (t„)\ : n = 0,1, ..,7V} < c\\u- x R \\ [h 1 + \\u - Xr\\ q ], 



(22) 



max ■ 



U n + ^ a iX[r 3 ,+oo)(tn) - X(t n ) 



= 0,1,..,7VV ^cWu-xrWo^ + Wu-XrWq]. (23) 
Proof : Similar to Theorem 4.6 of [3], based on Theorem 5.4 and Lemma 5.1. □ 



5.3 Superconvergence results for linear systems 

In this paragraph we suppose that f(x(t),t) = A(t)x(t) + b(t) where A{t) is a matrix of dimension 
dxd whose elements are measurable and bounded on [0, T], and b(t) and the columns of A(t) belongs 
toH K+1 (0,T;E). 



Lemma 5.2. Let w e H K+2 (I n ;E) such that w(t) + A T w(t) = 0,w(t~) = E n . The 



,W 



w 



(') 



O.ii 



<c\E n \ (l = 0,...,K+l), 
<ch 1/2 \E n \ (l = 0,...,K + 2). 



(24) 
(25) 



Moreover, let v G V K (I n \E) such that v(t) = E n — J. n Px(w)(r)dr, where Pk is the L 2 -projection 
onV K (I n :E), then 



IMIco,„<(l+^)|£n|, 

||«-Hl co .„<^+ 2 | J B„| 
< ch K+b ' 2 \E„ 



v — w\ 



v — w\\ 



0,n 



< ch K+3 / 2 \E n 



(26) 

(27) 
(28) 
(29) 
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Proof : See [3], Lemma 5.7. □ .„ 

Lemma 5.3. The following inequality is verified 

IlL^ + ^lklloJ < exp{cT^\\q\\ ). (30) 

Proof : See [3], Lemma 5.4. □ 

Theorem 5.6. Under the assumptions of Theorem 3.3 and if the impulses are state-independent, 
then for h small enough and u n € T jk+l (I n \ E) we have 

\U n -x R (t n )\< 1 + ch]/ 2 \\q\\ 0n |J7„_i-a;fl(t n _i)|+c||Q'|| 0|n ||u n -a;ii|| 0)n , (31) 

\\u n -x R \\ 0n < ch}J 2 \U n -i -x R (t n -i)\ +c\\u n -x R \\ Qn . (32) 

Proof : See [2], Theorem 4.10. D 

Remark 5.1. When the elements of A are in H (0,T; E), they are continuous and bounded, there- 
fore q is bounded and for any interval I n , \\q\\ n is proportional to h 1 ' 2 . D 



Theorem 5.7. Suppose that f(x(t),t) = A(t)x(t) + b(t) where b(t) and the columns of A(t) are in 



H +1 (0,T; E), and ^2 \ a j\ < +oo. Then if h is small enough we have 



je.J 

\U n - x R (t n )\ < (1 + ch) |f/ n _i - ccij(i„_i)| + ch K+3 / 2 \\u n - x R \\ , 

Proof : Similar to the proof of Theorem 5.10 of [2]. □ 

As a consequence of Theorems 5.1, 5.6 and 5.7 and Lemma 5.3 we obtain the next superconvergence 
result. 

Theorem 5.8. Under the assumptions of Theorem 5.7, we have the following result 
\U n -x R (t n )\<c\\U -x°\+h K+2 x^W V 7 Kl|. □ 

6 Numerical implementation 

In this paragraph we suppose that E = R, we shall also write V K (I n ) and V K (0, 1) for V K (I n \E) 
and V K ([0, 1]; E). We shall construct basis for P K (0, 1) and V K+1 (0, 1) and we deduce basis for 
V K (I n ) and T' K+1 (In), n = 1, .., N. Consider a (k + l)-points quadrature formula of the form 

,1 K+l 

/ ^(C)rfC= Va fc ^(%), (33) 

Jo fe=i 

where < r/i < ... < rji < rji + \ < ... < t]k+i < 1, such that this formula is exact for polynomials of 
degree at most IK + 1. Once the quadrature points have been obtained, we use them to construct 
bases for u n in P K (I n ) and v n in P K+1 {I n ). Let us denote by {<pk ■ k = 1, . . . , K + 1} the Lagrange 
interpolating polynomials associated with the K + l points {t;k}k=ii 

^(0=n-t 1 F-4' k = l,...,K + l, (34) 



and let 



4>nk{t) = 4>ki r^ j; tn-l<t<t n . (35) 
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.__ Then {4> nk ■ k = 1, . . . ,K + 1} is the desired basis of P K (I n ) 1 and the polynomial u n has the 
following representation 

EK+l 
U n k4>nk{t), (36) 

k—1 

where u nk = u n (t n k), t nk = tn-i + h n ^ k , k = 1, ..., if + 1. The basis for the polynomials i>„ 
in -P^ +1 (i„) is obtained in the following manner. We construct the new family of polynomials as 
follows 

Uk(T) = frMO<%, 0<r<l, k = l,...,K + l, 

\Vo(r) = l, 0<r<l. 

It is readily seen that the new family of polynomials {ip k : k = 0, . . . ,K+1} is a basis of P K+1 (0, 1). 
The set of polynomials 

i>nk{t)=il> k ( t ~* n ~ 1 ), t„-i<t<t„, (38) 

will be the desired basis for P K+1 (I n ). 

We now derive an equivalent set of equations for system (10). First, set v n (t) = ip n o(t) = 1 in (10). 
We obtain 



U n = U n -!+ / F(u n ,t)dt. (39) 

Then set v„ = ipnk, k = 1, ..., if + 1 in (10) and use the quadrature (33) to integrate exactly the 
polynomial u n ■ v n with u n given by (36) . We obtain 



U n tpnk(tn) I U n 1p nk (t)dt = U n -i1pnk(tn-l) + / F(u n ,t)lp nk (t)dt. 

1 JT+1 

But ip n k(t n ) = Vfe(l) = 0, and ^> n fc(* n _i) = Vfc(0) = J <l>k(O d C = E ai^kivi) - a k, then 



A"+] 

E 
1=1 

x+i 



/ u n (t)ijj nk (t)dt = ^ u nj - / 0„j(t) ( --^-ipnk(- — ^— -) ) dt, 



.3=1 
if+1 



1 ,; /* — i«-l 
/in 



K+l .1 

V u„j / <f>j{0<t>k{()d(, 

.7 = 1 J ° 



K+l K+l 



= - X! U «J X! O-j&jiV^&kiVi), 
j=l j=l 

V'nkQ'k- 



Also for fc = 1 , . . , if + 1 we have 

1 



Unk = U n - 1 + — F(u n ,t)ipnk(t)dt, (40) 



K+l 

with U n (t) = Y, Unk(t>nk{t)- 
k=\ 



The equation (10) is equivalent to: Uq = x°, and for n = 1, ..., N 



fc=i 
and 



K+i , , . . 
u n (t) = f/„_i + V ^^ / F(u n ,t)ip nk (t)dt 

zri a k Ji n 



U n = U n _ 1 + / F{u n ,t)dt, 
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where w n is the function defined by (9). 

If we assume the integrals containing the nonlinear term F are evaluated with the (K + l)-point 
quadrature (33) which is exact for polynomials of degree 2K + 1, then (10) leads to the following 
numerical scheme : Uq — x°, and 

h K+l 
u n k = U n -i + — y~] aiF(u n ,t n i)ipk(m)> 

a-k 

1=1 

for all k — 1 K + l and 
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j . .. j 



K+l 



U n = U n -i + h n y^ aiF(u n ,t n i). 



1=1 

K+l 

for all n = 1, ..., TV, with u„(t) = J^ u n k4>nk{t) and u„ is the function defined by (9). 

fc=i 

7 Numerical tests 

To validate the convergence results we consider two problems. In the first one, we consider a linear 
ordinary differential equation 

x(t) = 2tx(t), t e [0, 1] , 

x(0) = 1. 

while in the second one, we take a nonlinear ordinary differential equation 

x(t) = -2tx 2 (t), te [0,1], 
x(0) = 1. 

In this example the assumptions of Theorem 5.5 are satisfied with B = 2. In Figures 1 and 2 we 
consider the first example with one state-dependent impulse at the point 0.5 where ai(a;(0.5~)) = 
x(0.5~). Then we consider state-independent impulses. Two cases arc considered for each example. 
In the first case we suppose that we have 1000 impulses in t* = ~^r and a, = (— 1)* +1 . For the 
second case we suppose infinitely many state-independent impulses with an accumulation point at 
1/3, with 

1 1 , 1000 

Tt = — — -r—r and at = -—-—. 
3 3 l+1 2 1 

For each case, we compute the logarithm of L 2 -error (ln(eL2)) and the logarithm of the nodal error 
(ln(eN)) as functions of the logarithm of the step of the partition ln(h). 
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Figure 1: Nodal error with one state-dependent impulse 
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Figure 2: i 2 -error with one state-dependent impulse 
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Figure 3: Nodal error for example 1 with 1000 impulses 
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Figure 4: L 2 -error for example 1 with 1000 impulses 
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Figure 5: Nodal error for example 1 with infinitely many impulses 
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Figure 6: i 2 -error for example 1 with infinitely many impulses 
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Figure 7: Nodal error for example 2 with 1000 impulses 
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Figure 8: L 2 -error for example 2 with 1000 impulses 
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Figure 9: Nodal error for example 2 with infinitely many impulses 
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Figure 10: L 2 -error for example 2 with inhnitely many impulses 



Dubeau,Quansafi,Sakat 



REFERENCES 

111 P.G. Ciarlet, "The Finite Element Method for Elliptic Problems", North-Holland, Amsterdam, 
1978. 

M.C. Dclfour and F. Dubeau, "Discontinuous Polynomial Approximations in the Theory of One- 
Step, Hybrid and Multistep Methods for Nonlinear Ordinary Differential Equations", Math. 
Comp., 47 (1986), 169-189 and sl-s8. 

M.C. Delfour and F. Dubeau, "Fixed Mesh Approximation of Ordinary Differential Equations 
with Impulses", Numer. Math., 78 (1998), 377-401. 

[4] M. Delfour, W. Hager and F. Trochu, "Discontinuous Galerkin Methods for Ordinary Differential 
Equations", Math. Comp., 36 (1981), 455-473. 

F. Dubeau, "Approximation polynomiale par morceaux des equations differentielles" ', Ph.D. the- 
sis, Universite dc Montreal, Fevrier 1981. 

F. Dubeau, "Existence, Uniqueness and Approximation of the Solution of an ODE with (Infinitely 
Many) State-Dependent Impulses via a Fixed Mesh Galerkin Formulation" , Approx. Theory and 
its AppL, 15 (1999), 55-73. 

[7] F. Dubeau, "On First Order Ordinary Differential Equations with Infinitely Many State- 
Dependent Impulses", Differential Equations and Dynamical Systems, 5 (1997), 85-89. 

F. Dubeau, A. Ouansafi and A. Sakat, "Equations differentielles ordinaires avec impulsions qui 
dependent de l'etat", Annates des Sciences Mathematiques du Quebec, 24, no 2 (2001), 145-154. 

F. Dubeau, A. Ouansafi and A. Sakat, "Approximation dc la solution d'unc equation differcncicllc 
ordinaire avec impulsions qui dependent de l'etat" , CRM Technical Report no. 2706, Universite 
de Montreal, 2001. 

[10] J.L. Lions and E. Magenes, "Problemes aux limites non homogenes et applications", Vol. 1, 
Dunod, Paris, 1968. 

[11] A. Sakat "Existence, unicitc ct approximation de la solution d'une equation differencielle ordi- 
naire avec impulsions", Ph.D. thesis, Universite de Rabat, 2000. 



Approximation of State-Dependent... 



Journal of Computational Analysis and Applications,Vol. 7, No.4,437-446, 2005, Copyright 2005 Eudoxus Press, LLC 



The Newton Method for Hammerstein Equations 

J. A. Ezquerro and M. A. Hernandez 

University of La Rioja. Department of Mathematics and Computation. 

C/ Luis de Ulloa s/n. 26004 Logroho. Spain. 

E-mail: jezquer@dmc.unirioja.es, mahernan@dmc.unirioja.es 



Abstract 

We study the convergence of Newton's method to a solution of a nonlinear in- 
tegral equation of the Hammerstein type and second kind. We analyse the conver- 
gence of the method from a point of view different to traditional ones. We localize 
a solution of a particular integral equation and it is approximated directly by New- 
ton's method. Finally, some numerical tests are provided, where nonlinear integral 
equations of the Hammerstein type with degenerated kernels are considered. 

Keywords: nonlinear equations in Banach spaces, Newton's method, semilocal 
convergence theorem, recurrence relations, integral equation of the Hammerstein type. 
Classification A.M.S. 1991: 65J15, 47H17. 

The research reported herein was sponsored in part by the University of La Rioja 
(API-04/13) and MCyT (BFM 2002-00222). 

1 Introduction 

In this paper, we study a particular nonlinear integral equation of the Hammerstein type 
and second kind: 

x(s)=e(s)+ / !C(s,t,x(t))dt, (1) 

J a 

where K.(s,t,x(t)) = K(s,t)[x(t) 1+P + x(t) 2 ], p G [0,1], and K is continuous and non- 
negative in [a, b] x [a, b], and £ is a continuous function such that £(s) > 0, s G [a, b]. 

This type of integral equation has two interesting aspects for us. Firstly, if we use in 
the study of (1) the usual procedure, the method of successive approximations, we need 
that the kernel /C(s, t, x(t)) is Lipschitz continuous in its third argument (see [3]). Observe 
that (1) does not satisfy this condition. Moreover, in this paper, contrary to Davis, we 
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do not consider a homogenous situation. Secondly, it is known (see [4]) that by means 
of the study of Newton's method we can locate domains of existence and uniqueness of 
solutions for (1) if we consider: 

F:VC C[a, b] -+ C[a, b], V = {x e C[a, b]; x(s) > 0, s e [a, b}}, 

[F(x)](s) = x(s) - l(s) - [ K(s,t)[x(t) 1+p + x(t) 2 ]dt, pe[0,l], (2) 



and we set out F(x) = 0. 

The Newton method is the most used iteration to solve those equations as a conse- 
quence of computational efficiency even less speed of convergence can be got. Attention 
will now be given to the conditions under which the Newton sequence {x n } defined, start- 
ing from some point Xn S X, by the algorithm: 

x n+x =x n - F'(x n )~ 1 F(x n ), n > 0, (3) 

if (F'(a; n )) _1 = T n exists for n — 0, 1, 2, . . . , will converge to a solution x* of the equation 
F(x) = 0. 

Under the hypothesis of the Newton-Kantorovich theorem [6] : 

(Ci) ||r ||</3, 

(C 2 ) ||xi -Xo|| < 71, 



(C 3 ) ||F"(x)|| < L in some closed ball B(x ,p), 

(C 4 ) h = f3r]L< 1/2, 

the convergence of the Newton sequence {x n } implies the existence of the solution x*. This 
theorem gives not only conditions for the existence of x* but also information concerning 
the regions of existence and uniqueness of x* and error bounds for the terms x n of the 
Newton sequence as approximations to x*. 

Once the values of (3 and rj are calculated, we only need to determine ||F"(a;)|| in 
a neighborhood of xq. The original proof, given by Kantorovich [7], uses recurrence 
relations. In [6], a new proof is presented where the concept of majorant function is used. 

In [8], it is given another proof where majorant sequences are used and changes con- 
dition (C 3 ) by 

||F'(x)-i^)||<L||a;-3/||, L > 0, (4) 

where x, y are in an open convex domain Q C X. Other authors ([1], [2], [5]) consider the 
following generalization of the previous conditions: 

\\F'(x)-F'(y)\\<L\\x-y\\ p , L > 0, p6[0,l], (5) 
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where x, y are in an open convex domain OCX. Observe that if p — 1, (5) is reduced to 

(4). 

However, operator (2) satisfies neither (C3) nor (4) and (5) in T>. 
For (2), we have: 



[F\x)y](s)=y(s)- / K(s, t)[(l + p)x{tf + 2x(t)]y(t) dt, (6) 

J a 

[F"(x)yz](s) = - [ K(s,t)[(l+p)px(ty- 1 + 2]z(t)y(t)dt, 



and 

\\F"(x)\\ < A||s|| p - 1 + .B, xeV, (7) 

f b 

where A = (1 +p)pM, B = 2M and M = max / |i^(s,t)| dt for the max-norm. 

s6[o,6] J a 

The interest of this paper is to prove the convergence of the Newton method for 
operators with second Frechet derivative satysfying condition (7). For this, we use a new 
technique where a system of recurrence relations is constructed, and starting from them, 
the semilocal convergence of Newton's method is proved and a domain of existence of 
solutions of F(x) = 0, where F is given by (2), is also provided. After that, a domain of 
uniqueness of solutions for F(x) = is obtained. 

Finally, we consider a particular integral equation of type (1) and the domains of 
existence and uniqueness of solution are calculated. We also approximate its solution by 
the direct application of the Newton method. We finish with different numerical tests 
where integral equations of type (1) are considered. 

Throughout the paper we denote 



B(x, r) = {y G X; \\y — x\\ < r} y B(x, r) = {y G X; \\y — x\\ < r}. 

2 Localization of the solution 

Firstly, we locate a starting point xq for Newton's iteration, so that r = F'^o)" 1 exists. 
So, from (6), it follows that 

||/-F'(x)|| = ((l+p)\\x p \\ +2||x||)M. 

By the Banach lemma, if ((1 +p)||^ p || + 2||x||)M < 1, we have 

WF'ixV'W < 



I - ((1 + p)\\xP\\ +2||x||)M 

Therefore, r = -F^rrn) -1 G C(C[a, b], C[a, b}) exists for some xq G V, where C(C[a, b],C[a, b]) 
is the set of bounded linear operators from C[a, b] into C[a, b], if 

((l+p)||xg||+2||x ||)M <1. (8) 
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Moreover, 



|r °"- i-((i + p)KII + 2|M)m p - 



Furthermore, ||F(a;o)|| < \\xo — £\\ + (\\x p || + ||xg||) M, according to the definition of 
the operator F, and, in consequence, 

r w( M , ^ \\x - e\\ + (\\xi^\\ + \\xl\\) M 

\\Xl - X = r F(x ) < 77— ' ' p|| " ., 7TTTJ- = V- 

1- ((1+p) Hsgll +2||x ||)M 

The next aim is to prove the convergence of sequence (3) to a solution of the equation 
F(x) = 0, where the operator F is defined in (2), and to obtain the domains of existence 
and uniqueness of solution. For, the first, it suffices to see that (2) is a Cauchy sequence. 
Firstly, we introduce the following auxiliary equation in t: 

[2 - 3pt(A(\\x \\ - ty- 1 + B)]t - 2[1 - pt(A(\\x \\ - if' 1 + B)] V = 0, pe [0, 1], (9) 

and we assume that it has one positive root less than ||xo|| at least. The smallest root 
satisfying this condition is denoted by R (i. e.: R < ||xo||)- 
Secondly, for f(x) = (1 — a;) -1 , we observe that 

if a = (3R{A{\\x \\ - R) p - 1 + B) < 2/3. Therefore, it is supposed that a < 2/3. 
Thirdly, for x G B(xq, R), we have 

||/-r F'(x)|| < ||r ||||F , (x)-F , (xo)|| <«. 

Since a < 2/3 < 1, then, by the Banach lemma, F'(x)~ l exists and 

\\F\xy l F\x Q )\\<f{a). 

Moreover, for (10), \\xi — x \\ < rj < R, X\ G B(xq,R) and Fi = ^'(xi) -1 exists with 
||r 1 F'(:Fo)||</(a). 

From Taylor's formula and (3) it follows that 

r-x'i 



F( Xl ) = F(x ) + F'(x )(x 1 -x )+ / F"(x)(x 1 - x) dx 

JXQ 

F"(x + t(xi - x ))(x 1 - x )(l - t) dt (xi - x ) 2 , 



x 

1 

'0 

and consequently 

||^i)|| < / ||F"(xoH-i(xi-Xo))||(l-i)di||xi-x || 2 
Jo 
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< ^(A(\\x \\ - Rf- 1 + B)\\ Xl - xo|| 2 < ^(A(||x || - Rf- 1 + B)R\\ Xl - x \\. 
Thus 

Oi 

\W~x x \\ < ||riF'(xo)||||r ||||F(xi)|| <-j{a)\\x 1 -x^\\ < rj < R 

and 

||#2 - xo\\ < ||^2 - xi|| + ||xi — xo|| < f 1 + ^/(a)) rj < R. 

Hence we can follow the procdure because £2 £ B(xq,R). We then assume, for k = 
1, 2, . . . , n — 1, that the next conditions are true: 

(a) Tj, exists and ||rfeF'(xo)|| < /(«), 

(b) \\F(x k )\\< 1 -(A(\\x \\-Ry- l + B)R\\x k -x k ^\\, 

(c) ||a; fc+1 - x k \\ < -f(a)\\x k - x k ^\\, 

k 

(d) ||xfc+i-Xo|| < 77^(2 /(«)) < jR - 

i=0 

So, we can prove, by induction, that (a)-(d) also hold for k = n. 

As x n G B(xq,R), then r n exists and ||r„F'(a;o)|| < /(«)■ Since \\x n — x n -i\\ < R, we 
have 

||F(x n )|| < ^(A(\\x \\ - Rf' 1 + B)R\\x n - x n _i|| (11) 

and, by (10), we obtain 

||x n+ i-x n || < ||r n F , (xo)||||r ||||F(x n )|| < -/(a)||a: n -a: n _i|| < y^fi®)) \\xi-x \\ < R. 
Therefore 

n 

x \\ < \\x n+1 -x n \\ + ||x n -x || < V^2 (o 1 ^ V < R ' 

i=0 






71+1 



Finally, it follows immediately that (3) is a Cauchy sequence. Indeed, 

1 1 X n -\-m X n \\ \ ||X n _|_ m X n -\- rn — \\\ + ||X ra -|_ m _i X n _|_ m _2|| + ' ' " H~ H^n+l 2?n|| 

< L (2 /(a) J ll " 1 -"° 11 ^ l- f /(a) U /(a) J " 

and consequently {x n } is a Cauchy sequence, since a < 2/3. Thus, {x n } converges to a 
limit x* such that F(x*) = by letting n — ► 00 in (11). 

So, the folowing semilocal convergence theorem is given for the Newton method when 
it is applied to operator (2). 
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Theorem 2.1 We suppose that condition (8) is satisfied and let R be the small- 
est positive root of equation (9) which is less than \\x$\\. Assume that R exists, a = 
(3R(A(\\x \\ - R)p~ 1 + B) < 2/3, where p e [0, 1], 

A=(l+p)pM, B = 2M, j3 



l-((H-p)K||+2||x ||)M' 



and B(xq, R) C V. Then the Newton method (3) converges to a zero x* of operator (2). 



3 Uniqueness of the solution 

Now, we provide a result on the uniqueness of the zero x* of operator (2). 

Theorem 3.1 Under the hypotheses of the last theorem, we can guarantee the unique- 
ness of x* in B(xo,R). Moreover, x* is unique in V = B(x ,r) C\V, where r is the 
smallest positive root of the equation: 

(1 + p) P Af3\\x \\ p - 1 (t 2 - R 2 ) - 2Af3(t 1+p - R 1+p ) + (1 + p)p(t - R)(B{3(t + R) - 2) = 0, 

provided that r exists and r > R. 



Proof. Firstly, we prove that x* is unique in B(x 0} R). It is supposed that z* is 
another zero of (1) in B(xq,R). So, from the approximation 



= r o [F(z*)-F(x*)] 



J ToF'ix* +t(z* -x*))dt 
Jo 



[z — X 



P(z* - X* 



it follows that if the operator P = J Q T F'(x* + t(z* — x*)) dt is invertible, then z* = x*. 
Then, by the Banach lemma, we only have to prove taht \\I — P\\ < 1. Indeed, 

\\I-P\\ < ||r || / \\F'(x* + t(z*-x*))-F'(x )\\dt 
Jo 



<0 



i r i 



JO 



\F"(x + s(x* + t(z* - x*) - x ))\\ ds \\x* + t(z* - x*) - x \\ dt 



< (3(A(\\x \\ - Ry- 1 + B) / ((l-t)\\x*-x \\+t\\z*-x \\)dt 

Jo 

< P(A(\\x \\ - RY' 1 + B)R = a < 2/3 < 1. 
Secondly, we assume that z* is in T>$ = B(xq, r) fl T>. Similarly to the above, we have 



I/-PII < lir r 



/ F'(x* + t(z* -x*))- F'{x ) dt 
Jo 
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•1 r l 



F"{x + s((l - t)(x* - x ) + t(z* - x ))) ds ((1 - t)(x* - x ) + t(z* - x )) dt 

JO 

<A/3 / \\(x + s((l-t)(x* -x )+t(z* - xo))^ 1 ds {R + t{r - R)) dt 
Jo Jo 

•1 pi 



-B[3 f f {R + t{r-R))dsdt 
Jo Jo 



< Af3 f f (\\xo\\ p - 1 -sr- 1 \\(l-t)(x*-xo)+t(z*-x )\\ p - 1 ) ds (R + t(r - R)) dt 
Jo Jo 

<A(3 ! (R + t(r-R)) M|x || p - 1 --(i2 + t(r-i2)) p - 1 J dt + Bfi 1 -^ 
if r > R. ■ 



4 Application 

Next, we solve a nonlinear integral equation of type (1) by direct application of Newton's 
method. It is considered that K is a degenerated kernel, which is considered by other 
authors (see [3]), i. e., the kernel K is such that K(s,t) = S(s)T(t), so that the problem 
in infinite dimension can be solved. 

Let be the following particular integral equation of type (1): 

x(s) = 1 + s f t 8 [x(t) 3/2 + x(t) 2 } dt. (12) 



First of all, we calculate the domains of existence and uniqueness of solution. If we choose 
xq(s) = 1 for theorem 2.1, we have 

A = 1/12, 5 = 2/9, ^ = 18/11, 77 = 4/11. 

Equation (9) is then 

' ^-176+ ^548 + -JL=] t- ^264 + -2L ) t 2 ) = 



242 v V V^tJ V V^t 
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The smallest positve root and less than ||xo|| = 1 is R = 0.4289 .... So a = 0.2333 . . . 
and the hypotheses of theorem 2.1 hold. In consequence, (12) has a solution x* in 

{ueC[0,l\; \\u- 1|| < 0.4289...}. 

Following theorem 3.1, x* is unique in 

{u e C[0, 1]; ||u - 1|| < 5.3598 ...}DV. 

For (12), the operators F and F', given respectively by (2) and (6), are 

[F{x)]{s) = x(s) - 1 - s / t 8 [x{t) 3/2 + x{tf] dt, p G [0, 1], 

Jo 

[F'(x)y](s) = y(s) - s f ^[(S^)^^) 1 / 2 + 2x(t)]y(t) dt. 

Jo 

If it is supposed that (F'(x)) -1 exists, then y(s) = (F'(x))~ 1 w(s). So, we consider, 

y(s) = w(s) + si, 
where 

1=1 t s [(?,l2)x(t) 1 ' 2 + 2x(t)}y{t)dt. 
Jo 

If the last equality is multiplied by s 8 [(3/2)x(s) 1 / 2 + 2x(s)} and integrated between and 

1, we obtain 

_ jl s 8 [(3/2)a:(s) 1 / 2 + 2x(s)}w(s) ds 

~~ 1 - £ s 9 [(3/2)x(s) 1 /2 + 2x(s)]ds 

provided that L s 9 [(3/2)x(s) 1//2 + 2x(s)} ds ^ 1. Therefore 

M ,#, xx-i M M4 J 1 t 8 [(3/2) 3 ;(t)V 2 + 2 3 ;(t)] W (t)^ 

V (s = F (x)) wis) = wis) + s — = . 

l-^t^[(3/2)x(ty/ 2 + 2x(t)]dt 

The direct application of Newton's method is then 

x n +i(s) = x n (s) - (F'(x n ))~ 1 F(a;„)(s) = 1 + s -^— — ^- -, 

1 — ^"n 

where 

A n = [ t 8 [x n (t) 3/2 + x n (t) 2 ]dt, B n = [ t 8 [(3/2)x n (t) l/2 + 2x n (t)}x n (t)dt, 
Jo Jo 

C n = [ t 9 [(3/2)x n (t) 1/2 + 2x n (t)]dt, D n = [ t 8 [(3/2)x n (t) 1/2 + 2x n (t)]dt. 
Jo Jo 

We start at x (s) = 1 and the approximated solution is x*(s) = 1+0. 367498s (see table 1). 

Note that the approximated solution x* lies within the existence domain of solution 

obtained above (see figure 1). 
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% 


Xi(s) 


1 

2 
3 


1 + 0.341880s 
1 + 0.367558s 
1 + 0.367498s 



Table 1: Approximated solution of (12). 



1.4 



1.2 



x + R = 1.4289... 



x 



x = 1 



0.8 



0.6 



x - -R = 0.5711 .. . 



Figure 1: Graph of the approximated solution of (12) 
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